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Abstract 

The relativistic Vlasov-Maxwcll system describes the evolution of 
a collisionless plasma. The problem of linear instability of this system 
is considered in two physical settings: the so-called “one and one- 
half” dimensional case, and the three dimensional case with cylindrical 
symmetry. Sufficient conditions for instability are obtained in terms 
of the spectral properties of certain Schrodinger operators that act on 
the spatial variable alone (and not in full phase space). An important 
aspect of these conditions is that they do not require any boundedness 
assumptions on the domains, nor do they require monotonicity of the 
equilibrium. 
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1 Introduction 

We obtain new linear instability results for plasmas governed by the 
relativistic Vlasov-Maxwell (RVM) system of equations. The main un¬ 
knowns are two functions f ± = f±(t,x,v) > 0 measuring the density 
of positively and negatively charged particles that at time t € [0, oo) 
are located at the point x € and have momentum v £ R d . The 
densities / ± evolve according to the Vlasov equations 

or ± 

+v-S7 x f ± +F ± -V v f ± =0 (1.1) 

where v = v/\] 1 + |u| 2 is the relativistic velocity (the speed of light is 
taken to be 1 for simplicity) and where F^ 1 = F^t, x , v) is the Lorentz 
force, given by 


F ± = ± (E + B ext + vx{B + B ext )) 

with E = E(t, £c) and B = B (t,x) being the electric and magnetic 
fields, respectively, and E ext ,B e:ct possible fixed external fields. The 
self-consistent fields are retrieved from Maxwell’s equations 

<9B <9E 

V • E = p, V • B = 0, VxE = --, VxB=j + —, 

where 

P = p(t,x) = J(f+ - f~) dv (1.2) 

is the charge density and 

j=j {t,x) = J v{f + - f~) dv (1.3) 

is the current density. In addition to the speed of light, we have taken 
all other constants that typically appear in these equations (such as 
the particle masses) to be 1 so as to keep the notation simple. 

Novelty of the results. Let us mention the main novel aspects 
of our instability results: 

Unbounded domains: our problems are set in unbounded domains 
(as opposed to domains with boundaries or periodic domains). One 
consequence is that the spectrum of the Laplacian (which shall appear 
prominently) has an essential part. 

Non-monotone equilibrium: we do not assume that the equilibrium 
in question is (strongly) monotone (see (11.711 below). Many estimates 
in previous works rely heavily on monotonicity assumptions. 
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Construction of equilibria: in Isection 7l we construct nontrivial 
equilibria in the unbounded, compactly supported 1.5 d case. Previously, 
this has been done in the periodic setting by means of a perturbation 
argument about the trivial solution which is a center (in the dynamical 
systems sense). The proof here relies on a fixed point argument. 

1.1 Main results 

For the convenience of the reader, we provide the full statements of our 
results here, although some necessary definitions are too cumbersome. 
We shall refer to the later sections for these definitions. 


The physical setting. As is explained in detail below we consider 
two problems: the 1.5 dimensional case and the 3 dimensional case 
with cylindrical symmetry. We shall refer to these two settings as the 
1.5d case and the 3d case , respectively, for brevity. In a nutshell, we 
consider these settings because they provide enough structure so that 
basic existence and uniqueness results hold and because they possess 
well-known conserved quantities which may be written explicitly. 


The equilibrium. The conserved quantities mentioned above - 
the microscopic energy e ± and momentum - are the subject of 
further discussion below (see (11.171) for the 1.5d case and (11.281) for 
the 3d case), however we stress the fact that they are functions that 
satisfy the time-independent Vlasov equations. Hence any functions of 
the form 

f°’ ± (x,v)=p, ± (e ± ,p ± ) (1.4) 

are equilibria of the corresponding Vlasov equations (for positively 
or negatively charged particles). The converse statement - that any 
equilibrium may be written in this form is called Jeans’ theorem [T. 
In Isection 71 we show how to obtain an actual equilibrium of the full 
Vlasov-Maxwell system in both physical settings. When there is no 
room for confusion we simply write /r ± (e,p) or pA instead of p ± (e ± ,p ± ). 
We shall always assume that 

0 < / 0,± (a;, v) e C 1 have compact support fl in the ^-variable. 

Again, the existence of such equilibria is the subject of IsectionTl In 
addition, we must assume that 

there exist weight functions uA = c(l + |e ± |) _ “ ^ 

where a > dimension of momentum space, and c > 0, 


such that the integrability condition 


5p± 

de 



(e ± ,p ± )<w ± (e ± ) 


( 1 . 6 ) 


4 


















holds. This implies that J (\p.f\ + |/ip |) dxdv < oo in both the 1.5d 
and 3 d cases, where we have abbreviated the writing of the partial 
derivatives of pA 1 . This abbreviated notation shall be used throughout 
the paper. It is typically [2] assumed that 

fif < 0 whenever pA > 0. (1.7) 

We call this a strong monotonicity condition. We do not make any 
such assumption. 

The main results. To facilitate the understanding of our main 
results we state them now, trying not obscure the big picture with 
technical details. Hence we attempt to only extract those aspects of the 
statements that are crucial for understanding, while refering to later 
sections for some additional definitions. First we define our precise 
notion of instability: 

Definition 1.1 (Spectral instability). We say that a given equilibrium 
pA is spectrally unstable, if the system linearised around it has a purely 
growing mode solution of the form 

(e At / ± ( a b «), e At E(a:), e At B(x)) , A > 0. (1.8) 

We also need the following definition: 

Definition 1.2. Given a self-adjoint operator (bounded or unbounded) 
A, we denote by neg(_4) the number of negative eigenvalues (counting 
multiplicity) that it has whenever there is a finite number of such eigen¬ 
values. 

In our first result we obtain a sufficient condition for spectral in¬ 
stability of equilibria in the 1.5d case. The condition is expressed in 
terms of spectral properties of certain operators that act on functions 
of the spatial variable alone. 

Theorem 1.1 (Spectral instability: 1.5dcase). Let / 0,± (x, v) = p ± (e,p) 
be an equilibrium of the 1.5 d system (11.151) satisfying (II.Gil . There exist 
self-adjoint Schrodinger operators A 5 and A 2 and a bounded operator 
B° (all defined in (11.241 ) ) acting only on functions of the spatial variable 
(and not the momentum variable) such that the equilibrium is spectrally 
unstable if 0 is not in the spectrum of A ? and 

neg (A 0 2 + ( B 0 )* {A^y 1 B° ) > neg (A 0 ,) . (1.9) 

The second result provides a similar statement in the 3 d case with 
cylindrical symmetry, as discussed in further detail in Isubsection 1.41 
below. 
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Theorem 1.2 (Spectral instability: 3d case). Let / 0 ,± (a;, v) = /x ± (e,p) 
be a cylindrically symmetric equilibrium of the RVM system satisfying 
«ni). There exist self-adjoint operators A^A® an d ^3 an d a bounded 
operator B( (all defined in (11.331 ) ) acting in the spatial variable alone 
(and not the momentum variable) such that the equilibrium is spec¬ 
trally unstable if 0 is not an L e -eigenvalue of A® (see 
below), 0 is not an eigenvalue of A® (0 will always lie in the essential 
spectrum of Ai, but this is not the same as 0 being an eigenvalue) and 

neg + (£°) * (I?) B?) > neg (l?) + neg (l°) . (1.10) 

Let us make precise the notion of an L 6 -eigenvalue. 

Definition 1.3 (L 6 -eigenvalue). We say that X GS. is an L 6 -eigenvalue 
of a self-adjoint Schrodinger operator A. : C L 2 (R";R m ) — > 

Z, 2 (R";R m ) given by A = —A + K, if there exists a function 0 ^ 
u x G H( oc (S. n ;R m ) n L 6 {R n ;R m ), with Vn A G L 2 (R n ;R m ) n , such 
that Au x = A u x in the sense of distributions. The function u x is 
called an L°-eigenfunction. 

Remark 1.1. We remark that L e is a natural space to consider in 
three dimensions due to the compact embedding H 1 ^) L 6 (Q) where 

XI C R 3 is a bounded and smooth domain. In fact, any function which 
decays at infinity and whose first derivatives are square integrable, also 
belongs to L 6 (R 3 ). Therefore this is a natural condition for the po¬ 
tential formulation of Maxwell’s equations where there is no physical 
reason for the potentials to be square integrable, but the condition that 
the fields are square integrable corresponds to the physical condition 
that the electromagnetic fields have finite energy. 

The proofs of these two theorems appear in Isubsection 4.11 and 
Isubsection 4.21 receptively. Let us describe the main ideas of the 
proofs. For brevity, we omit the ± signs distinguishing between posi¬ 
tively and negatively charged particles in this paragraph. Since we are 
interested in linear instability, we linearise the Vlasov equation around 
f°. The only nonlinear term is the forcing term F ■ V„/, so that the 
linearisation of m becomes 

+ w ■ V x / + F° • V„/ = — F ■ V„/° (1.11) 

at 

where F () is the equilibrium self consistent Lorentz force and F is the 
linearized Lorentz force. We make the following growing-mode ansatz: 


Definition 1.3 


Ansatz: the perturbations (/, E, B) have 

time dependence e xt , where A > 0. 


( 1 . 12 ) 


Equation (11.111) can therefore be written as 


(X + V) f = —F • V„/° 


(1.13) 


6 











where 


X) — v • Vjc + F° • 


(1.14) 


is the linearised Vlasov transport operator. We then invert (11. Kill by 
applying A (A + V) 1 . which is an ergodic averaging operator along the 
trajectories of V (depending upon A as a parameter), see Q] Eq. (2.10)]. 
Hence we obtain an expression of / in terms of a certain average of 
the right hand side —F • V„/° depending upon the parameter A (see 
(E2D and G2D)- This expression for / is substituted into Maxwell’s 
equations through the particle and charge densities, resulting in a sys¬ 
tem of (elliptic) equations for the spatial variable alone (recall that the 
momentum variable is integrated in the expressions for p and j). The 
number of linearly independent equations is less than one would expect 
due to the imposed symmetries. However, in both cases the equations 
can be written so that they form a self-adjoint system denoted A4 a 
( see (13.511 for the 1.5d case and (13.1411 for the 3d case) that has the 
general form 


M X 


—A + 1 

0 


0 

A- 1 


+ JC X 


acting on the electric and magnetic potentials, where 7C A is a uniformly 
bounded and symmetric family. 

The problem then reduces to showing that the equation A\ x u = 0 
has a nontrivial solution for some value of A > 0. The difficulty here 
is twofold: first, the spectrum of J\4. x is unbounded (not even semi- 
bounded) and includes essential spectrum extending to both +oo and 
—oo. Second, for each A, the operator A4 a has a different spectrum: 
one must analyse a family of spectra that depends upon the parameter 
A. In [3] we address the following related problem: 


Problem 1.1. Consider the family of self-adjoint unbounded operators 


M X = A+JC X 


—A + 1 

0 

+ 

7C A + 

ic\_- 

0 

A - 1 

k A + 

/C A _ 


A G [0,1] 


acting in L 2 (R d ) ® L 2 (IR d ), where {7C A }Ae[o,l] a uniformly bounded, 
symmetric and strongly continuous family. Is it possible to construct 
explicit finite-dimensional symmetric approximations of A4 a whose 
spectrum in (—1,1) converges to that of M. x for all A simultaneously? 

A solution to this problem allows us to construct finite-dimensional 
approximations to Ai X . We discuss this problem in Isubsection 2.21 
The conditions (11.91) and (11.101) appearing in the main theorems above 
translate into analogous conditions on the approximations, and those, 
in turn, guarantee the existence of a nontrivial approximate solution. 
Since the approximate problems converge (in an appropriate sense) to 
the original problem, this is enough to complete the proof. A crucial 
ingredient is the self-adjointness of all operators: this guarantees that 
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the spectrum is restricted to the real line. The strategy is to “track” 
eigenvalues as a function of the parameter A and conclude that they 
cross through 0 for some value A > 0. To do so, we require knowledge 
of the spectrum of the operator A1 a for small positive A, which is 
obtained from the assumptions (USD and (11.101) . and for large A which 
arises naturally from the form of the problem. 

Yet even with a solution to this problem at hand, some difficul¬ 
ties remain. In the cylindrically symmetric case there is a geometric 
difficulty. Namely, cylindrical symmetries must be respected, a fact 
that requires a somewhat more cumbersome functional setup. In par¬ 
ticular, the singular nature of the coordinate chart along the axis of 
symmetry requires special attention. To circumvent this issue we shall 
do all computation in Cartesian coordinates, and use carefully chosen 
subspaces to decompose the magnetic potential. The second difficulty 
is the lack of a spectral gap, which is due to the unbounded nature of 
the problem in physical space. As a consequence, the dependence of 
the spectrum of Ai x on A is delicate, especially as A —> 0, and needs 
careful consideration. 

1.2 Previous results 

Existence theory. The main difficulty in attaining existence re¬ 
sults for Vlasov systems is in controlling particle acceleration due to the 
nonlinear forcing term. Hence existence and uniqueness has only been 
proved under various symmetry assumptions. In [5] global existence 
in the 1.5d case was established and in @] the cylindrically symmetric 
case was considered. Local existence and uniqueness is due to El- 

Stability theory. One of the important early results on (linear) 
stability of plasmas is that of Penrose [TJj. Two notable later results 
are IP. We refer to [T] for additional references. The current result 
continues a program initiated by Lin and Strauss mm and continued 
by the first author mm- In |21 m the equilibria were always as¬ 
sumed to be strongly monotone, in the sense of (H3- This added sign 
condition (which is widely used within the physics community, and is 
believed to be crucial for stability results) allowed them to obtain in [9| 
a linear stability criterion which was complemented by a linear instabil¬ 
ity criterion in m- Combined, these two results produced a necessary 
and sufficient criterion for stability in the following sense: there ex¬ 
ists a Schrodinger operator C° acting only in the spatial variable, such 
that C° > 0 implies spectral stability, and C° ^ 0 implies spectral 
instability. In mm the monotonicity assumption was removed, which 
mainly impacted the ability to obtain stability results. The instability 
results are similar to the ones of Lin and Strauss, though the author 
only considers the 1.5d case with periodicity. This is due to his meth- 
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ods which crucially require a Poincare inequality. We remark that our 
results recover all previous results when one restricts to the monotone 
case. 

1.3 The 1.5 d case 

First we treat the so-called 1.5d case, which is the lowest dimensional 
setting that permits nontrivial electromagnetic fields. In this setting, 
the plasma is assumed to have certain symmetries in phase-space that 
render the distribution function to be a function of only one spatial 
variable x and two momentum variables v = (iq,^), with iq being 
aligned with x. The only non-trivial components of the fields are the 
first two components of the electric field and the third component of 
the magnetic field: E = (Ei, E 2 , 0) and B = (0,0, B), and similarly for 
the equilibrium fields. The RVM system becomes the following system 
of scalar equations 


ft/ ± + fiia x / ± ± {Ei + d 2 b)^ 1 / ± ± {e 2 

-v 1 B)d V2 f ± = 0 (1.15a) 

dtEi = —ji 

(1.15b) 

d t E 2 + d x B = j'2 

(1.15c) 

o. 

II 

(1.15d) 

dtB = -d x E 2 

(1.15e) 


where p and ji, J 2 are defined by ED and ED- 

1.3.1 Equilibrium 

I 11 [section 71 we prove the existence of equilibria f°' ± {x, v), which can 
moreover be written as functions of the energy and momentum 

f°’ ± (x,v)=p, ± {e ± ,p ± ) ( 1 . 16 ) 

as in ED, where the energy and momentum are defined as: 

e± = (v) ± (j>° (x) ± (j) ext (x), p ± = V 2 ± V’°( x ) ± ip ext { x ) ( 1 - 17 ) 

and where (v) = -^/l + |t;| 2 , and <p° and 4’° are the equilibrium electric 
and magnetic potentials (both scalar, in this case), respectively: 

d x ct>° = -E°, d x 4>° = B° ( 1 . 18 ) 

and similarly (f> ext and 4> ext are external electric and magnetic poten¬ 
tials that give rise to external fields Ef xt and B ext . It is a straightfor¬ 
ward calculation to verify that and are conserved quantities of 
the Vlasov flow, i.e. that V±e ± = 'D±p ± = 0, where the operators T>± 
are defined below, (11.201) . 
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1.3.2 Linearisation 


Let us discuss the linearisation of (11.151) about a steady-state solution 
(/°' ± ? E°,B 0 ). Using ansatz (11.121) and Jeans’ theorem (11.41) . together 
with (11.171) and (11.181) the linearised system becomes: 

(A + V±)f ± = Ei ± nf viB =f {pfv 2 + Pp)E 2 (1.19a) 


XEi = ~j i (L19b) 

A E 2 + d x B = -j 2 (1.19c) 

d x E 1 = p (1.19d) 

A B = -d x E 2 , (1.19e) 

where 


V ± = v x d x ± + Et xt + v 2 (B° + B ext ))d Vl ± (E° - i, 2 (B° + B ext ))d V2 

(1.20) 

are the linearised transport operators as in (11.141) . and 

P = J (f + - f~)dv, ji = J Vi{f + - f~) dv (1.21) 

are the charge and current densities, respectively. 

We now construct electric and magnetic potentials (j) and tp, re¬ 
spectively, as in (11.181) . Equation (ll,19bl) implies that E\ has the same 

spatial support as j\ which is a moment of which, in turn, has the 

same x support as ^ (this can be seen from (I1.19al) for instance). We 
deduce that E\ is compactly supported in fl and choose an electric 
potential </> € H 2 (Q) such that Ei = —d x (j) in fl and Ei = 0 outside 
Q. Since E\ vanishes at the boundary of Q, we must impose Neumann 
boundary conditions on </>, and as Ei depends only on the derivative 
of cj) we may impose that (f> has mean zero. The magnetic potential ip 
is chosen to satisfy B = d x ip and E 2 = —A ip (this is due to (I1.19cl) l. 
Then the remaining Maxwell’s equations (I1.19bl) - (|1.19dl) become 

A d x <p = —XEi = ji in Q (1.22a) 

(—dx + \ 2 )ip = —d x B — \E 2 =j 2 in R (1.22b) 

— d^.(p = d x Ei = p in Q (1.22c) 

where (ll.22cl) is complemented by the Neumann boundary condition 

—d x (p = Ei = 0 on dCl. 

The linearised Vlasov equations can now be written as 

(A + V±)f ± = ±nfv! d x (p ± PpVid x ip ± A {pfv 2 + Pp)ip 2g 
= ±pfV±(p ± pfV±ip ± A (nfv 2 + Pp)ip 

where we have used the fact that T>±u = V\d x u if u is a function of x 
only. 
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Now let us specify the functional spaces that we shall use. For the 
scalar potential <j> we define the space 

L&Sl) := {/G L 2 (n) : jf / = o| 

while for the magnetic potential we simply use L 2 (K), the stan¬ 
dard space of square integrable functions. We denote by H k ( R) (resp. 
H k (fil)) the usual Sobolev space of functions whose first k derivatives 
are in L 2 (R) (resp. L 2 (fi)). Moreover, we naturally define 

H k (n) := {/G H k (n) : ^/ = oj 

and the corresponding version which imposes Neumann boundary con¬ 
ditions 

H k rl (n.) := { / G H k (Q) : d x f = 0 on Oft} . 

Finally, to allow us to consider functions that do not decay at infin¬ 
ity we use the conditions (LSD and m to define weighted spaces 
£± as follows: we take the closure of the smooth and compactly sup¬ 
ported functions of (x, v) (with the x support contained in 12) under 
the weighted-L 2 norm given by 

||u||c + = [ w ± \u\ 2 dvdx 

J nxR 2 

and we denote the inner product by (•, -) £ . In particular we can view 
any function u{x) £ L 2 (VL) or Lq(SY) as being in £± by considering u 
as a function of (x, v) which does not depend on v. We can extend 
this to functions in L 2 (R) by multiplying them by the characteristic 
function of the set f2, denoted In- Hence the function In itself may 
be regarded as an element in £±. 

1.3.3 The operators 

Finally, we define the operators used in the statement of IThcorcm 1.11 
First define the following projection operators: 

Definition 1.4 (Projection operators). We define Q± to be the or¬ 
thogonal projection operators in C± onto ker(2?±). 

Remark 1.2. Although this definition makes reference to the spaces 
£±, the operators Q± do not depend on the exact choice of weight 
functions w±. This may be seen by writing {Q\h){x,v) as the point- 
wise limit of ergodic averages along trajectories (see \Remark S.l\ and 
\Lemma 6.1\ ) . 

This allows us to define the following operators acting on functions 
of the spatial variable x, not the full phase-space variables: 

A°h = — d 2 h + f Y, — 1 )hdv (1.24a) 

J ± 
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A 2 h = —dlh— j [ipV 2 dv^j h - J ^t) 2 M±Q± [«2 h]dv 


B°h= ( dv ) h+ / i^h] dv 


(B°)*h= (j £>±dt^/i + J J2fitv 2 Q±hdv. 


(1.24b) 

(1.24c) 


(1.24d) 


Their precise properties are discussed in Isubscction 6.11 For future 
reference, we mention the important identity 



(1.25) 


which is due to the fact that = yfv 2 


dp- 


1.4 The cylindrically symmetric case 

Since notation can be confusing when multiple coordinate systems are 
in use, we start this section by making clear what our conventions are. 


Vector transformations and notational conventions. We 

let x = (x, y, z) = xei + ye 2 + ze 3 denote the representation of the 
point x £ R 3 in terms of the standard Cartesian coordinates. We 
define the usual cylindrical coordinates as 

r = \J x 1 + y 2 , 9 = arctan (y/x), z = z. 

and the local cylindrical coordinates as 

e r =r~ 1 {x,y, 0), e g =r~ 1 (-y,x, 0), e z = (0,0,1). 

By cylindrically symmetric we mean that in what follows no quantity 
depends upon 9 (which does not imply that the 9 component is zero!). 
When writing f(x) we mean the value of the function / at the point 
x in Cartesian coordinates. We shall often abuse notation and write 
f(r, 9, z) to mean the value of / at the point (r, 9 , z) in cylindrical coor¬ 
dinates. A point v £ R 3 in momentum space shall either be expressed 
in Cartesian coordinates as 

v xyz = (v x ,v v , v~) = (v ■ ei)ei + (v ■ e 2 )e 2 + (v ■ e 3 )e 3 

or in cylindrical coordinates (depending upon the point x £ R 3 in phys¬ 
ical space) as 

VrzB = (v r , v e ,v z ) = (v ■ e r )e r + (v ■ e g )e B + (v ■ e z )e z . 
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However we shall not be too pedantic about this notation, and shall 
use v (rather than v xyz or v rz g) when there’s no reason for confusion. 

A vector-valued function F shall be understood to be represented in 
Cartesian coordinates. That is, unless otherwise said, F = ( F x , F y ,F z ) = 
F x e i +F y e 2 +F z e 3 . Its expression in cylindrical coordinates shall typ¬ 
ically be written as F = F r e r + Fgeg + F z e z . 

Differential operators. Partial derivatives in Cartesian coordi¬ 
nates are written as d x , d y and d z , while in cylindrical coordinates they 
are d r ,dg and d z . They transform in the standard manner. It is im¬ 
portant to note that since we work in phase space, we shall require 
derivatives with respect to v as well. One important factor appearing 
in the Vlasov equation is v x ■ V x , which transforms as 

(v x ■ Va,)/i = v x d x h -(- Vydyh H- v z d z h 

= v r d r h + r~ 1 vgdgh + v z d z h 
= v r d r h + r^ 1 vg{vgd Vr h — v r d vg h) + v z d z h. 

However the next term in the Vlasov equation transforms “neatly”: 

(F • V„)/i = F x d Vx h + F y d Vy h + F z d Vz h 

= ( F x cos 6 + F y sm6)d Vr h + {—F x sind + F y cos 6)d vg h + F z d Vz h 
= F r d Vr h + Fgd Vg h + F z d Vz h. 

1.4.1 The Lorenz gauge 

As opposed to the system (11.191) , here we do not get a system of scalar 
equations. We begin with Maxwell’s equations. It is well known that 
there is some freedom in defining the electromagnetic potentials <p (we 
use ip in the cylindrically symmetric case rather than (f) to avoid con¬ 
fusion) and A, satisfying 

d t A + \7ip = -~E, V X A = B. 

Remark 1.3. Whenever the differential operator V appears without 
any subscript, it is understood to be V x , that is, the operator (d x ,d y , d z ) 
acting on functions of the spatial variable in Cartesian coordinates. 

The same holds for the Laplacian A. 

We choose to impose the Lorenz gauge V ■ A + ^ = 0, hence 
transforming Maxwell’s equations into the hyperbolic system 

— ip-Atp = p, (1.26a) 

o 2 

— A-AA=j. (1.26b) 

We remark that this is not unique to the cylindrically symmetric case, 
and the expressions above are written in Cartesian coordinates. 
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1.4.2 Equilibrium and the linearised system 

We define the steady-state potentials ip° : R 3 — > R and A 0 : R 3 — > M 3 
as 

V</?° = — E°, VxA° = B° 

which become 

E° = —d r ip°e r — d z ip°e z , B° = — d z A° g e r H— d r (rAg)e z . (1.27) 

r 

The energy and momentum may be defined (analogously to (11.171) 1 as 
e% l = (v)±tp°(r,z)±<p ext (r,z), 

(1.28) 

p% l =r(v e ±A° g (r,z)±A? t (r,z)), 

where we remind that ( v ) = \/l + |'y xy , z | 2 . It is straightforward to 
verify that they are indeed conserved along the Vlasov flow (which is 
given by the integral curves of the differential operators T >±, defined in 
(11.301) below). To maintain simple notation we won’t insist on writing 
the cyl subscript where it is clear which energy and momentum are 
meant. The external fields are also assumed to be cylindrically sym¬ 
metric and their potentials satisfy equations analogous to (11.271) . We 
recall OD, namely that any equilibrium is assumed to be of the form 


f°’ ± (x,v)=v ± (e ± ,P ± )- 


Considering the Lorenz gauge, and applying the ansatz (11.121) and 
Jeans’ theorem m the linearisation of the RVM system about a 
steady-state solution (/^.E^B 0 ) is 

(A + V±)f ± = ±(A + V±)(nf<p + r/j,p (A • e e )) ± A nf(-<p + A • v) 


(1.29a) 

\ 2 tp — A<p = j(f + - /“) dv (1.29b) 

A 2 A — AA = J (/ + — f~)v dv (1.29c) 

where 

V± = v xyz ■ V x ± (E° + E“* + v xyz x (B° + B ext )) ■ V v 

= v r d r + V z d z + (±E° r ± E e r xt ± Vg(B° z + B ext ) + W 1 ^^)^ 

+ (±v z (B° r + B e r xt ) T v r (B° z + B ext ) + r~ 1 vgv r )d ve 
± (. E° z + E e z xt + v e (B° r + B* xt ))d Vz 

(1.30) 

are the linearised transport operators. The Lorenz gauge condition 
under the growing mode ansatz is 


V • A + A (p = 0. 


(1.31) 
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1.4.3 Functional spaces 

Even more so than in the 1.5 d case, choosing convenient functional 
spaces is crucial, due to the singular nature of the correspondence 
between Cartesian and cylindrical coordinates. We define 

L%i (® 3 ) = the smallest closed subspace of L 2 (R 3 ) com¬ 
prised of functions which have cylindrical symmetry. 

A short computation using cylindrical coordinates shows that the de¬ 
composition L 2 (R 3 ) = L 2 cyl { R 3 ) © (L^^R 3 ))- 1 - reduces the Laplacian. 

This means that the Laplacian commutes with the orthogonal pro¬ 
jection of L 2 (R 3 ) onto L 2 yl ( R 3 ). Hence the Laplacian is decomposed 
as 

A = A cy i + A cyl ± • 

As we have no use for (L^^R 3 )) 1 - we shall abuse notation slightly and 
denote A cy i as simply A. We now consider vector valued functions 

Ae^,(K 3 ;K 3 ):= (i^(R 3 )) 3 - 

We decompose such functions as 

A = (A ■ e g )e g + ((A ■ e r )e r + (A ■ e 2 )e z ) 

= A e + A rz e L^(R 3 ;R 3 )0L 2 ,(R 3 ;R 3 ). 

By computing with cylindrical coordinates we once again discover that 
this reduces the vector Laplacian A on L 2 yi (R 3 ;R 3 ). Note that this 
reduction does not occur for A on L 2 (R 3 ;R 3 ) (i.e. without the cylin¬ 
drical symmetry). 

We further define the corresponding Sobolev spaces H^ yl (R 3 ),Hg (R 3 ; R 3 ), 
iL^(M 3 ;R 3 ) of functions whose first k weak derivatives lie in L 2 yl ( R 3 ), 
Lg(R 3 ;R 3 ) and L 2 ,(R 3 ;R 3 ), respectively. Note that, because of the 
reductions above, A is self-adjoint on L 2 yl ( R 3 ) with domain H 2 yl { R 3 ) 
and A is self-adjoint on each of Lg(R 3 ;R 3 ) and L 2 Z (R 3 ;R 3 ) with do¬ 
mains Hg( R 3 ;R 3 ) and H 2 Z (R 3 ;R 3 ), respectively. 

As in the 1.54 case we shall require certain weighted spaces that 
allow us to include functions that do not decay. We define as the 
closure of the smooth compactly supported functions m : R 3 x R 3 -> R 
which are cylindrically symmetric in the x variable, and have a;-support 
contained in f2, under the norms 

ll u lk+ = / w ± \u\ 2 dvdx, 

where the weight functions u ^ are the ones introduced in m- 

1.4.4 The operators 

We now define the operators used in the statement of IThcorcm 1.21 As 
in the 1.54 case we shall require the following definition of projection 
operators: 
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Definition 1.5 (Projection operators). We define Q± to be the or¬ 
thogonal projection operators in onto ker (T>±). 

As in the 1.5d case, the operators do not depend upon the ex¬ 
act choice of weights . Now we are ready to define the operators 
of the cylindrically symmetric case. For brevity, given v = (v r ,vg,v z ), 
we define vg = vgeg and v rz = v r e r + v z e z . All operators act on 
functions of the spatial variables only: the operator A? acts on func¬ 
tions in L^, yl ( ffi 3 ), A ° on functions in Lg(R 3 ;K 3 ), A 3 on functions in 
L^(K 3 ;R 3 ), and Bi on functions in Lg(R 3 ;K 3 ) with range L^ yZ (R 3 ). 
We have: 

Aih = —Ah+ f ^ — l)hdv (1.33a) 

J ± 


A° 2 h = -A h ~ y J dvj h - J '^vep±Q°±\h ■ v e \ dv 



(1.33b) 

A°h=-Ah- 1 ^ Vrznt Q± [h ' Wz] dv 

J ± 

(1.33c) 

Bi h = f ^p.t(.Q\-l)[h-v e ]dv 

J ± 

(1.33d) 

(B?)Vi = [ J2l l tvs{Q° ± -l)hdv 

(1.33e) 

The precise properties of these operators are discussed 
We also mention an identity analogous to (11.251) 

inlsubsection 6.21 

J (rpf + vg/lf)dv = 0 

(1.34) 


which is due to the integrand being a perfect derivative: ^ = r/z± + 
vgpf. 


1.5 Organization of the paper 

fn lsection 2l we provide some necessary background, including the cru¬ 
cial result on approximating spectra found in [3]. Then we treat the two 
problems - the 1.5d and 3d cases - in parallel: in lsection 31 we formu¬ 
late the two problems as an equivalent family of self-adjoint problems 
which we then solve in Isection 4l The proofs of the main theorems 
are concluded in Isection 51 In Isection fil we provide the rigorous treat¬ 
ment of the various operators appearing throughout the paper, and in 
Isection 7l we show that there exist nontrivial equilibria. 


16 


















2 Background, Definitions and Notation 


In this section we remind the reader of the various notions of con¬ 
vergence in Hilbert spaces in order to avoid confusion. For a Hilbert 
space ft we denote its norm and inner product by ||-||^ and (•,•)%, 
respectively. When there is no ambiguity we drop the subscript. We 
denote the set of bounded linear operators from a Hilbert space ft to a 
Hilbert space © as *8 (fj, ©), and when ft = © we simply write (fj). 
The operator norm is denoted |H| where, again, when there is no 
ambiguity we may drop the subscript. 

Definition 2.1 (Convergence in *8(ij, ©)). Let T, T n € 93 (fj,©), 
where ngN. 

(a) We say that the sequence T n converges to T in norm (or uni¬ 
formly) as n —► oo iff \\T n — —> 0 as n —> oo. In this case 

we write T n —t T. 

(b) We say that the sequence T n converges to T strongly as n —> oo 
iff we have the pointwise convergence T n u —> Tu in ft for all 
u £ ft. In this case we write T n A T. 

Now let us recall some important notions related to unbounded 
self-adjoint operators: 

Definition 2.2 (Convergence of unbounded operators). Let A and A n 
be self-adjoint , where n £ N. 

(a) We say that the sequence A n converges to A in the norm resolvent 
sense as n —>■ oo iff (A n — z)^ 1 —> (A — z ) _1 for any z£C\I. 
In this case we write A n "’A A. 

(b) We say that the sequence A n converges to A in the strong re¬ 
solvent sense as n —»■ oo iff (, A n — z)~ x A (A — z)~ 1 for any 
zeC\R. In this case we write A n A. 

Remark 2.1. Notice that for any self-adjoint operator A, the resolvent 
(A — z)^ 1 is a bounded operator for any z £ C \ R. 

2.1 Basic facts 

The subsequent results will be used throughout the paper without 
explicit reference. 

Lemma 2.1. Let ft, © be Banach spaces, T,T n £ 93 (if, ©) and u, u n £ 
where n £ N, and assume that T n A T and u n —> u. Then T n u n — t u 
as n —> oo. 

Proof. We compute 


|| T n U n - Tu\\ < || T n {u n - w)|| + ||(7A — T)u\\ 
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This supremum is finite by the uniform boundedness principle, so the 
first term converges to zero since The second term converges 

to zero since T n T. □ 

Corollary 2.1. If T n T and S n A S as n —> oo then T n S n —> TS 
as n —> oo. 

The following result complements Weyl’s theorem (see [5} IV, The¬ 
orem 5.35]) on the stability of the essential spectrum under a relatively 
compact perturbation. In our setting we know more about the pertur¬ 
bation than being merely relatively compact. 

Lemma 2.2. Let A = —A + K. : H 2 (R n ) C L 2 (R") ->• L 2 (R ra ) be 
a self-adjoint Schrddinger operator with K. £ ©(L 2 (ffi. ra )) and 1C = 
KIP where V : L 2 (R n ) —> L 2 (R") is the multiplication operator by the 
characteristic function 1q of some bounded domain fl C R". Then A 
has a finite number of negative eigenvalues. 

Proof. We want to show that the set 

A f := {u : ( Au,u ) < 0} = : ( K.u,u ) < — ||Vu|| 2 | 

is finite dimensional. Consider the larger set 

NcN':={u: ||/C|| ||Pu|| 2 > ||Vu|| 2 }. 

We claim that it is sufficient to restrict to functions u supported in 
Cl. Indeed, letting \ be a non-negative, compactly supported smooth 
cutoff function which equals 1 in Cl, we have 

<-Au,u> = ||Vu|| 2 > ||xVu|| 2 = ||[Xj V]u + V(xu)|| 2 

= II V(x«)H 2 + II [X, V]u || 2 + 2 <[x, V]«, V(X«)) 

> ||V(xu)|| 2 -2||Vx|| L oo ||u|| ||V(x«)|| 

>^l|V(x«)|| 2 -C||u|| 2 

where [x, V]/ = \S7f — V(x/) = (Vx)/ is the commutator between 
V and multiplication by x- By choosing x appropriately, the constant 
C may be made arbitrarily small. This proves the claim. Finally, by 
the boundedness of Cl the set J\P is clearly finite dimensional, as there 
can only be finitely many energy levels below any given threshold in a 
bounded domain. □ 

2.2 Approximating strongly continuous families of 
unbounded operators 

Here we summarise the main results of [3j on properties of approxi¬ 
mations of strongly continuous families of unbounded operators. We 
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shall require these results in the sequel. We refer to [3] for full details, 
including proofs. Let = ,f) + ©fj_ be a (separable) Hilbert space with 
inner product (•,•} and norm ||-|| and let 


A x 


A\ 0 

0 -.4 A 


and /C A 


X\ + K, X J 

IC X + JC X _ 


A £ [0,1] 


be two families of operators on depending upon the parameter A £ 
[0,1] (the range [0,1] of values of the parameter is, of course, arbi¬ 
trary), where the family _4 A is also assumed to be defined for A in an 
open neighbourhood Do of [0,1] in the complex plane. They satisfy: 


i) Sectoriality: The sesquilinear forms corresponding to 
are sectorial for A £ Do, symmetric for real A, have dense domains 
£)(a±) independent of A £ D 0 , and fl 0 3 Ai-> a±[it, u] are holomorphic 
for any u,v £ 2)(a^.). [In the terminology of 0, a^. are holomorphic 
families of type (a) and A x are holomorphic families of type (B).\ 


ii) Gap: A± > 1 for every A £ [0,1]. 

iii) Bounded perturbation: {/C a }a£[o,i] C (ij) is a symmetric 
strongly continuous family. 


iv) Compactness: There exist symmetric operators V± £ 18 (fj±) 
which are relatively compact with respect to the forms a^., satisfying 
K, x = KAV for all A £ [0,1] where 


V = 


V+ 

0 


0 

V- 


Finally, if the family A x does not have a compact resolvent we as¬ 
sume: 


v) Compactiflcation of the resolvent: There exist holomorphic 
forms {in^jAeDo of type (a) and associated operators {W±}ag.Do of 
type (B) such that for A £ [0,1], W A are self-adjoint and non-negative, 
and if tu A is the form associated with 


W A 


'W A 0 

0 -W A ’ 


A £ Dq, 


then D(tu A )n£>(a±) are dense for all A £ D 0 and the inclusion (J)(m A )n 
0(a), ||• IIa A ) (-ft) INI) is compact for some A £ D 0 and all e > 0, 

where a A is the form associated with 

A x :=A x + eW x , A £ Dq, e > 0. (2.1) 
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Define the family of (unbounded) operators {Ad A }Ae[o,i]i acting in 
h, as 

M X = A X + 1C X , A e [0,1]. 

Our main result in [3] is: 

Theorem 2.1. Let A x be as in m, and define 
M X =A X +JC X , A e [0,1]. 

Let ( e efc}fceN C Sj be a complete orthonormal set of eigenfunctions 
of A x , let Q x n : Sj -4 Sj be the orthogonal projection operators onto 
span(e A 1; ..., e x n ) and let A4 X n be the n-dimensional operator defined 
as the restriction of M x to Q x n (Sj). Fixe* > 0, and define the function 

E : [0,1] x [0, £*] —> ( subsets of (— 1,1), dn) 

E(A, £ ) = (-l,l)nsp(M A ) 
and for fixed e > 0 the function 

E e : [0,1] x M 4 ( subsets of (—1,1), du) 

E e (A,n) = (-l,l)nsp(M A „) 

where dn is the Hausdorff distance, N = N U {oo} and where we use 
the convention that Ai x ^ := Ai x . Then E and E e are continuous 
functions of their arguments in the standard topologies on K (and its 
subsets) and N. 

We recall the definition of the Hausdorff distance between two 
bounded sets S,TcC: 

dfi{ H, T) = max sup inf \x — y\, sup inf \x — y\ 

VyeT^GH xesv^ r 

3 An equivalent problem 

3.1 The 1.5 d case 

We will now reduce the linearised Vlasov-Maxwell system (11.191) to a 
self-adjoint problem in i 2 (R) x Lg(fi) depending continuously (in the 
norm resolvent sense) on the parameter A > 0. 

3.1.1 Inverting the linearised Vlasov equation 

Rearranging the terms in (11.231) we obtain 

(A + V±)f ± = ±(A + V±)(pf(f> + n±ijj) ± A/z^(-</> + v 2 ip), (3.1) 

where we use the fact that are constant along trajectories of the 
vector-fields T>±. In order to obtain an expression for f ± in terms of 
the potentials (f>, ip we invert the operator (A + T>± ) and to do this we 
must study the operators D±. 
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Lemma 3.1. The operators D± on £± satisfy: 

(a) T>± are skew-adjoint and the resolvents (A + T>±)~ 1 are bounded 
linear operators for Re A 0 with norm bounded by 1/| ReA|. 

(b) T>± flip parity with respect to the variable v\, i.e. if h(x, v±, v^) € 
T)('D±) is an even function of v\ thenT>±h is an odd function of 
v± and vice versa. 

(c) For real A 0 the resolvents ofT>± split as follows: 

(A + T>±)~ 1 = A(A 2 - Vl)- 1 - V±(X 2 - Vl)- 1 

where the first part is symmetric and preserves parity with respect 
to v\, and the second part is skew-symmetric and inverts parity 
with respect to v\. 

Proof. Skew-adjointness follows from integration by parts, noting that 
are in the kernels of T>± (to be fully precise, only skew-symmetry 
follows. However, skew-adjointness is a simple extension, see e.g. m 
and in particular exercise 28 therein). The existence of bounded re¬ 
solvents follows. The statement regarding parity follows directly from 
the formulas for T>± term by term. Finally, for the last part we use 
functional calculus formalism to compute 

1 \-V ± A V± 

A + V± ~ A 2 - V 2 ± ~ A 2 - X>| _ A 2 - V 2 ± ' 

As T>± are skew-adjoint, T>± are self-adjoint and hence the first term 
is self-adjoint and the second skew-adjoint. For the parity properties 
we note that as D± flip parity, T>± preserve parity and hence so do 
A 2 — T>± and their inverses. □ 

Applying (A + V±)~ l to (13.11) yields, 

= ±/i^±/i^±A(A + 2?±)’ 1 [^(-0 + i)2i/’)]- (3.2) 

Furthermore, using iLennna 3. 11 we split f± into even and odd functions 
of v\: 


ft, = ±pt<l > ± ± M?A 2 (A 2 - D±) 1 [-</> + V2lf] 

ftd = T^A2?±(A 2 - Vlr^-cf + vrf] 

using the fact that (j), if and p are all even functions of v\. For brevity, 
we define operators Q± : £± —> £± as 

Q± = A 2 (A 2 - £>±)“\ A > 0. 

When A —> 0 their strong limits exist, and are defined in IDcfinition 1.41 
(this convergence is proved in ILennna 6.11) . 


21 












Remark 3.1. Operators Q± also appeared in the prior works 
0/. In each of these Q± were defined as an integrated average over the 
characteristics of the operators D±. In fact, as the Laplace transform 
of a semigroup is the resolvent of its generator we see that the operators 
Q± in these prior works have the rule: 

/ O r-oo 

Xe Xs e sV± hds = A / e- Xs e~ sV± hds = A(A + V±)~ 1 h. 

-oo J 0 

Here we have defined the operators Q± directly from the resolvents 
of D± as this makes some of its properties clearer, although both ap¬ 
proaches have advantages. In particular we are able to split X{T>± + 
A)^ 1 into symmetric and skew-symmetric parts in \Lemma 3.1\ which 
simplifies some computations. 


3.1.2 Reformulating Maxwell’s equations 

Now we substitute the expressions ( 13 . 21 ) into Maxwell’s equations ( 11 . 221 ) . 
This shall result in an equivalent system of equations for </> and ip. Due 
to the integration dv we notice that f° d and f± d v 2 both integrate to 
zero, so that p and j'2 only depend on . 

Remark 3.2. It is important to note that due to the continuity equa¬ 
tion it is possible to express either (ll.22al) or (ll.22bl) using the remain¬ 
ing two equations in ( 11 . 221 ) . See \Lemma 5.4\ 

Gauss’ equation (ll.22cl) . Gauss’equation becomes 

-9x4> = J(fev ~ fev) dv 

(vt# + Mp + Q±[f^(-4> + V 2 lp)]) dv 

± 

+pt^)4’dv+ I ^utiQi-l^-cp + fa^dv, 

± J ± 

(3.3) 

where we have pulled pf outside the application of Q± as they belong 
to ker(2?±). 
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Ampere’s equation (jl.22b|) . Similarly, Ampere’s equation be¬ 
comes 


(~d 2 x + = J V 2 Utv - fev) dv 

= / Xb 2 + Vp + Q±[Mt(-<!> + V 2^)]) dv 
■' ± 

= / +tjiv 2 )4)dv 

+ J E^2 pj{Q± - l)[-(t>+V24>\dv. 


(3.4) 


An equivalent formulation. We write the two new expressions 
(13.31) and (13.41) abstractly in the compact form 


¥ 


+ \ 2 4> - n 


'o' 

¥ 


d 2 x (f) + p 


0 


(3.5) 


where, for A > 0, A4 a is a self-adjoint matrix of operators mapping 
L 2 (K.) x L\ j(f2) —> L 2 (R) x L§(fl) fsee lLcnmia 6.41) . We claim that this 
operator may be written either as 


M x = 


-a 2 + A 2 O' 
0 dl 


- J 


or, equivalently, as 


M x = 


a x (z x y 

B x -A$ 


(3.6) 


(3.7) 


where the various operators appearing above are given by 


J X 


h 

. 9 . 



± 1 + v\ 
(v) 




+E 


v 2 

-1 


v£(Q±~ !) 




(3.8a) 


A x h = -d 2 h + J J2nf(Q± - l)hdv 


(3.8b) 


A x h = -d 2 h + \ 2 h - 


\ ± 


J HpV 2 dv^j h- y^E^2 pf 2± I®2 h) dv 


B x h = 


iyj J2dpdv^jh + J E^2± [v 2 h\dv 


(3.8c) 


(3.8d) 
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(3.8e) 


(B x )'h 



nfv2Q±hdv. 


Remark 3.3. Though A > 0 in the foregoing discussion, all operators 
can he defined for A = 0, as we have already done for some (see (11.241 ) ). 

The expression (13.71) is no more than a rewriting of (13.31) and (13.41) . 
However the expression (13.61) requires some attention. In particular, to 
obtain it one has to use (11.251) as well as the integration by parts 


/ 


dv 2 


V2 dv 



± 1 + v'[ 
(v) 


dv. 


The properties of the operators appearing in (13.81) are discussed in 
details in lLcmma 6.21 and ILcmma 6.31 Let us briefly summarise: 

• Ai : C Ll(n) I/q(Q) is self-adjoint and has a purely 

discrete spectrum with finitely many negative eigenvalues. 

• Aff ■ H 2 (M.) C T 2 (R) —> L 2 (R) is self-adjoint, has essential spec¬ 
trum in [A 2 , oo) and finitely many negative eigenvalues. 

• B x : L 2 (R) —> Lq(Q) is a bounded operator, with bound inde¬ 
pendent of A. 

• : L 2 (R) x Lp(12) —> L 2 (R) x Lq(Q) is a bounded symmetric 
operator, with bound independent of A. 


3.2 The cylindrically symmetric case 

Our approach here is fully analogous to the one presented in lsubsection 3.11 
hence we shall keep it brief, omitting repetitions as much as possible. 

For convenience we denote analogous operators by the same letter, but 
we shall add a tilde to any such operator in this section. Hence, e.g. 
the operators analogous to T>± shall be denoted T>±. 


3.2.1 Inverting the linearised Vlasov equation 

Recall the linearised Vlasov equation (11.29al) 

(A + V±)f ± = ±(A + T)±)(nf(p + rp±( A ■ e g )) ± A + A • v). 
Inverting, we get the expression 

f ± = ± r/Xp (A ■ e e ) ± pf A(A + V ± )~ l {-ip + A • v), (3.9) 

and, recalling that we only care about the quantity / + — /“, we write 
it for future reference: 

/+-/- + A-e e ) + ^/x±A(A + 22 ± ) _1 (-^ + A-i)). 

± ± ± 

(3.10) 
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Lemma 3.2. The operators D± on satisfy: 

(a) T>± are skew-adjoint and the resolvents (A + T>±)~ 1 are bounded 
linear operators for Re A 0 with norm bounded by 1/| ReA|. 

(b) V± flip parity with respect to the pair of variables ( v r ,v z ), i.e. if 
h £ T)(V±) is an even function of the pair ( v r ,v z ) then V±h is 
an odd function of(v r ,v z ) and vice versa (see \Remark 3.4\ below). 

(c) For real A 0 the resolvents ofT>± split as follows: 

(A + V ± )~ l = A(A 2 - Vl)- 1 - V±{X 2 - V\Y l (3.11) 

where the first part is symmetric and preserves parity with respect 
to ( v r ,v z ), and the second part is skew-symmetric and inverts 
parity with respect to (v r ,v z ). 

We leave the proof, which is analogous to the proof of lLcmma 3.11 
to the reader. 

Remark 3.4. For a function h expressed in cylindrical coordinates as 
h(x,v r ,v z ,vg), we say that h is an even function of the pair ( v r ,v z ) 
if h(x, v r , v z , vg) = h(x,—v r ,—v z ,vg), where we flip the sign of both 
variables simultaneously. Note that this is a weaker property than both 
being an even function of v r and an even function ofv z . Odd functions 
of (v r ,v z ) are defined similarly. 

As in the 1.5d case, we define averaging operators. However, in this 
case both the symmetric and skew-symmetric parts are required. The 
operators Q± iSyrn and Q± skew map 9T± to 9t± and are defined by the 
rules 

Q±,sym = A 2 (A 2 — A > 0 

Ql, skew = -ND ± Y 2 A > 0. 

Note that by (13.111) we have A(A + V±)~ 1 = Q± ySym + Q± iSkew - 

3.2.2 Reformulating Maxwell’s equations 

We now rewrite Maxwell’s equations (ll.29b[) - (|1.29cD as an equivalent 
self-adjoint problem using the expression (13.101) . We start with (ll.29bl) : 

0 = A V -Aip- j (/+ - f~) dv 

= \ 2 ip - Ay - f (/^V> + rpp (A ■ eg) + pf\(X + V ± )~ 1 (-^ + A ■ v 

(3.12) 
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where ip £ Sj v . Next, the system of equations (11.29cl) becomes 


0 = A 2 A - AA - y (/+ - f~)v dv 

= A 2 A - AA - / ^ (vtv + r Bp ( A ■ e e) + M±A(A + V ± )~ 1 {-p + A • v)J vdv, 

(3.13) 

where A = (Ag,A rz ) £ Lg(R 3 ; R 3 ) x L 2 Z (R 3 ;M 3 ). As in (13.51) . we shall 
write these equations as a single system of the form 


A@ 


'O' 

<P 

= 

0 

A-rz. 


0. 


(3.14) 


which is a self-adjoint operator in io, see lLcmma 6.81 In analogy with 
(ED, we define 


A 


M 


A 

(sir 

(B^r 

B$ 


-(*$)■ 



-As 


(3.15) 


With v = (v r ,vg,v z ), we recall the notation vg = vgeg and v rz = 

—A 

v r e r + v z e z introduced before. Then the components of At are now 
given by 


A X 1 h=-Ah+\ 2 h+ 1 J2vt(Q±,sym- 
J ± 

- 1 )hdv (3.16a) 

A%h = — A h + X 2 h — J ^ PpVg dv^j 

h - / d 0ptQ±, S ymi h ' «e] dv 
■' ± 

(3.16b) 

A^h = -Ah + X 2 h - / ^ v rz pf Q± :Sym [h ■ v rz ] dv (3.16c) 

■' ± 

B x h= 1 J2Bt(Q X ±,sym-m-vg\dv 
J ± 

(3.16d) 

(B x )*h = 1 J2BtMQ±, sym -l)hdv 
■' ± 

(3.16e) 

&2 h = J ^2d't i} rzQ±.skewi h ' Vg\dv 

(3.16f) 


± 
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(B X )*h = - [ J2vtteQ±,ske W [Vrz • h\ dv 

J ± 

(3.16g) 

&3 h = j ^vfvrzQi'Skew h dv 

(3.16h) 


± 


{BzYh 



VrzQ 


A 

±,skew 


[h ■ v rz ] dv. 


(3.16i) 


These are derived from (13.121) and (13.131) . where some terms vanish 
due to parity in (v r ,v z ), (see ILcnmia 3.2f cW In particular, in every 
occurrence of A(A + V±)~ 1 = Q± ySym + Q± yS kew > exactly one of these 
operators vanishes after integration dv. In addition, we have made use 
of (11.341) . We further define an operator J as 


~ A 

J = 


A 2 -A 
0 
0 


0 0 

-A 2 + A 0 

0 -A 2 + A 


-At 


Let us briefly discuss these operators in further detail (their precise 
properties are treated in lsubsection 6.2|) : 

• The operators 

^ : H 2 cyl ( R 3 ) C Ll yl (R 3 ) -»■ L 2 cyl (R 3 ) 

A% : 7/f (R 3 ;R 3 ) C Lj{ R 3 ;R 3 ) L^(R 3 ;R 3 ) 

A% : i7 r 2 2 (M 3 ;R 3 ) C L 2 2 (R 3 ;R 3 ) -> L 2 2 (R 3 ;R 3 ) 

are self-adjoint, have essential spectrum in [A 2 ,oo) and a finite 
number of eigenvalues in (—oo, A 2 ). 

• The operators 

B^(R 3 ;R 3 )->L^(R 3 ) 

: L 2 0 ( R 3 ;R 3 ) L 2 2 (R 3 ;R 3 ) 

&3 : B 2 y i(M. 3 ) —»■ Z/ 2 .(R 3 ;R 3 ) 
are bounded, with bound independent of A. 

• J X :L 2 (R 3 ;R 3 )xL 2 yZ (R 3 )xL 2 2 (R 3 ;R 3 ) ^L 2 (R 3 ;R 3 )xL 2 yi (R 3 )x 
L 2 2 (R 3 ;R 3 ) is a bounded symmetric operator with bound inde¬ 
pendent of A. 


4 Solving the equivalent problem 

The problem is now reduced to finding some A £ (0,oo) for which 
the operators Ai x (in the 1.5 d case) and Ad (in the cylindrically 
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symmetric case) have non-trivial kernels (not the same A in both cases, 
of course). Our method is to compare their spectrum for A = 0 and 
A very large, and use spectral continuity arguments to deduce that as 
A varies an eigenvalue must cross through 0 (we shall show that both 
operators are self-adjoint, see ILcmma. 6.4l and lLcmma 6.8l below. hence 
the spectrum lies on the real axis). 

4.1 The 1.5 d case 

4.1.1 Continuity of the spectrum at A = 0 

Recall the condition (11.91) which we require for instability: 

neg(Al° + (23°)* K)" 1 !? 0 ) > negK). (4.1) 

We wish to move this condition to values of A greater than 0: 

Lemma 4.1. Assume that m holds and that zero is in the resolvent 
set of A®. Then there exists A* > 0 such that for all A € [0, A*] 

neg(Al 2 A + (S a )*(AI a )- 1 0 a ) > neg(Al A ). 

Proof. The proof follows immediately from the following three simple 
steps: 

Step 1. *4 a is invertible for small A > 0. We know from 
ILcmma 6.31 (below) that „4 A is continuous in the norm resolvent sense 
and has discrete spectrum. The norm resolvent continuity implies that 
its spectrum varies continuously in A, so as 0 is not in its spectrum 
at A = 0 there exists A* such that 0 is not in the spectrum for 0 < 
A < A*. Hence for all such A, Af is invertible and the operator „4 2 + 
(B x )*(Ai)~ 1 B x is well defined. 

Step 2. neg(„4. A ) = neg(.4.5) f° r a h A € [0, A*]. The spectrum 
of A x is purely discrete and 0 is in its resolvent set. This means that 
none of its eigenvalues can cross 0 for small values of A. 

Step 3. neg(A x + {B x )* (A^ 1 B x ) > neg(^ + (23°)* pf?)" 1 # 0 ) 
for all A € [0, A*]. Observe that 

• [0, oo) 3 A i-3 Aff + (B X )*(A X )~ 1 B X is norm resolvent continuous, 

• Al 2 + (B x )*(A x )~ l B x has essential spectrum in [A 2 ,oo), 

• Al 2 + (B x )*(A x )- 1 B x has finitely many negative eigenvalues. 

These statements follow from arguments similar to those appearing in 
the proof of ILemma 6.3l aT(c~). the last by the boundedness of the per¬ 
turbation and the location of the essential spectrum (see ILcmma 2.21) . 
Since 0 is not in the resolvent set at A = 0 we pick a < 0 larger than all 
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the (finitely many) negative eigenvalues of + (B 0 )* (Ai)~ 1 B° . The 
continuous dependence of the spectrum (as a set) on the parameter 
A implies that for small values of A no eigenvalues cross a and the 
number of negative eigenvalues can only grow as A increases. □ 


4.1.2 Truncation 


We follow the plan hinted at in Isubsection 2.21 first we discretise the 
spectrum, then truncate. The only continuous part in the spectrum of 
M x is due to Aff , hence we let W ( x) be a smooth positive potential 
function satisfying W(x) —> 00 as x —> ±00 which we shall add to A%. 
It is well known that the Schrodinger operator — <9 2 + W on L 2 (R) is 
self-adjoint (on an appropriate domain therein) with compact resolvent 
(and therefore discrete spectrum). Moreover, Cq°( R) is a core for 
both df + W and d 2 . Thus our approximating operator family is 
}a£[a*,oo),££[ 0 ,oc)> where 


M x = 


-A x e (B x r 


'-<9 2 + eW 

O' 

B x -A$ 


0 

<9 2 


A 2 O' 

0 0 


-J x 


■M 


defined on L 2 (R) x Lq(SY) and where A* is as given in lLcmma 4.11 For 
e > 0 this operator has discrete spectrum. As indicated in the state¬ 
ment of IThcorcm 2.1l we define truncated versions using the eigenspaces 
of the operator A. x . As this operator is diagonal, we can choose the 
eigenvectors to lie in exactly one of L 2 (R) or Lq(12). We denote the 
nth truncation, a projection onto an eigenspace of dimension 2 n con¬ 
sisting of n eigenvectors in each of L 2 ( R) and Lq(12), as M. X n which 
is self-adjoint and defined for e > 0, A > 0, n £ N. Moreover, the 
mapping A sp(A4^„) is continuous (that is, the set of eigenval¬ 
ues varies continuously). In particular if there are A* < A* for which 
neg(A4g* n ) ^ neg(A4g n ) then there must exist \ e , n € (A*, A*) for 
which 0 € sp(A4g n ). We have therefore just proved: 

Lemma 4.2. Fixe > 0,n € N. Suppose that there exist 0 < A» < A* < 
00 such that neg(A4 e * n ) ^ neg(A4 E „). Then there is a A Ei „ € (A*, A*) 
for which ker(A4 Era ) is non-trivial. 

The next step is thus to establish estimates on neg(A i X n ). 


4.1.3 The spectrum for large A 

We begin by looking at neg(Ai x n ) when A is large. This turns out to 
be relatively simple due to the block form of the untruncated operator. 

Lemma 4.3. There is A* > 0 such that for all A > A*, e > 0, and 
n £ N, the truncated operator Ai x n has spectrum composed of exactly 
n positive and n negative eigenvalues. In particular neg(A4 e n ) = n. 
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Proof. Take u = (ui,0) <£ L 2 (R) xLg(fl) with m G S(^ iE ), ||«|| L 2 (R)xL 2 (n) 
1 and m in the 2 n dimensional subspace associated with the truncation. 
Then, 



(”4l,e,ra' u l; u l)£2( R ) (-4l,e u lj u l )£2( R ) 


Wi) i2 ( K ) +£ 


V^UI 


2 

L 2 (R) ‘ 


As the second term is non-negative we may apply ILcnmia 6.3l df to 
see that, for all large enough A (independently of n and e), A4 A n 
is positive definite on a subspace of dimension n, so has n positive 
eigenvalues. Performing the same computation onti= (0,«2) in the 
subspace associated with the truncation and with u 2 G S(-A A £ ), we 
obtain that for large enough A, A4 a n is negative definite on a subspace 
of dimension n. As AA. Z n has exactly 2n eigenvalues the proof is 
complete. □ 


4.1.4 The spectrum for small A 

We now consider sp(A4 A * ( ). We recall the result on spectra of real 
block matrix operators in [2]: 

Lemma 4.4. Let M be the real symmetric block matrix 


M = 


A 2 B t 
B —A\ 


with Ai invertible. Then M has the same number of negative eigen¬ 
values as the matrix 


N = 


A 2 + B T Af 1 B 0 

0 —Ai 


Lemma 4.5. Assume that (14.11) holds and that zero is in the resolvent 
set of Af Then there exist A*,e* > 0 such that for all e G (0, £*) there 
is N > 0 such that for all n > N the operator A4 A * n satisfies 

neg(A4 A *„) > neg(A° + (B°)*(.A?) _1 S 0 ) + n - neg(A°). 

Proof. The number A* is the one given in ILcmma 4.11 and satisfies 
that for all A G [0, A*] the kernel of _4 A is trivial. Since eigenvalues 
(counting multiplicity) are stable under strong resolvent perturbations 
(see [H VIII.3.5.Thm 3.15.]), there exists e* > 0 such that neg(A 2 i£ + 
(S°)*(Ai) _1 S 0 ) > neg(^1° + OB 0 )*^?)- 1 # 0 ) for all e e [0,e,]. The 
result then follows from ILernma 4.41 since neelAt A l) = n eel A 2*e + 
(S A *)*(^*)" 1 B A *) + n-neg(yf A *). □ 
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4.2 The cylindrically symmetric case 


For brevity we write 


A 


M 


^ (*£)*' 

~A ~ A 

IS 4 — *<4.4 


where 


'A (Bp* 
B$ 


and B 4 


BJ 

B\ 


4.2.1 Continuity of the spectrum at A = 0 

Lemma 4.6. Assume that (11.101) holds, that A 3 does not have 0 as an 
L e -eigenvalue (see \Definition I73\ ) and that A 4 does not have 0 as an 
eigenvalue. Then there exists A* > 0 such that for A € [0, A*], 

neg(^ + {B^)* pP^Bp > neg(.A 4 ). 

Proof. We first note that as the mean perturbed charge is zero, it fol¬ 
lows from direct computation on the Green’s function of the Laplacian 
(see [9] Lemma 3.2]) that any L-eigenfunction of A® will also be square 

integrable and so be a proper eigenfunction. Note also that £>3 = 0. 

~o - 

Thus our assumptions imply that A 4 has no L -eigenfunction of 0. 


We model the proof on that of ILemma 4.11 splitting it into 4 steps. 


Step 1. A 4 is invertible for small A > 0 when restricted to 

functions supported in Cl. Let ^L 2 yi (R 3 ) x L 2 _(R 3 ; R 3 )^ 

be multiplication by the indicator function of Cl. We claim that for all 

small enough A > 0, V{A 4 )~ 1 V is a well defined bounded operator 

that is strongly continuous in A > 0 and has a strong limit as A —► 0. 

To prove this, we argue that if this were not the case, then zero would 

be an L 6 -eigenvalue, a contradiction. 

As L^ y/ (M 3 ) x L 2 ,(R 3 ;R 3 ) is a closed subspace of L 2 (R 3 ;R 4 ) we 

may work in the larger space to ease notation. To this end, let ||-|| and 

(•,•) denote the L 2 (R 3 ;R 4 ) norm and inner product. We can express 
~A 

A 4 in the form 

A 4 u = — A u + A 2 u + K x u 

where 7C A is uniformly bounded, strongly continuous in A > 0 and 

7C A = VK X V. 


Step 1.1. A 4 is bounded below when restricted to 
functions supported in Cl. First we claim that there exist con¬ 
stants A' > 0 and C > 0 such that we have the uniform lower bound 


lnAl 4 tt A 


> G||lnM A 


VA g (0, A'] 


(4.2) 
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where the constant C does not depend on A or on u x and where u x 
~ a , 

satisfies *4 4 « = 0 outside fl. Indeed, if not there would be sequences 
A ra —> 0 , it„ with ||loMn || L 2 = 1 that satisfy 

*4-4 u n = ~A .11 + A + /C^ 77 u n = f n —> 0 (4-3) 

as n —> 00 , with f n supported in Q. Hence, 

\\ U nf + I|Vm„|| 2 + ( K. Xn U n , U n ^ = ( f n ,U n ) 0 

so that || Vm „|| 2 is uniformly bounded for large enough n. Hence there 
exists a subsequence (we abuse notation and keep the same sequence) 
such that Vu„ v weakly in L 2 (R 3 ;R 4 ) for some v £ L 2 (R 3 ;R 4 ). 
By the standard Sobolev inequality ||<^|| i e(R3) < C || V<^|| L2 ( R 3) we have 
u n —^ u weakly in L 6 (R 3 ;R 4 ) for some u £ L 6 (R 3 ;R 4 ). Furthermore, 
by Rcllich’s theorem we have u n —> u in L 2 oc (R 3 ;R 4 ). This implies 
that necessarily v = S7u. In particular we deduce that the limit u is 
in L 6 (R 3 ;R 4 ) and ||lfi«|| = 1 so « ^ 0. Passing to the limit in ( 14 . 31 ) . 
u satisfies 

—Am + /C°m = 0 

in the sense of distributions and by elliptic regularity u £ Hf a C (R 3 ; R 4 ). 
In fact u is an L 6 -eigenfunction of 0, which we assumed impossible. 
This proves the claim. 


Step 1.2. A. v is invertible for all small enough A > 0. 

~ A 

For any A > 0, 0 does not lie in the essential spectrum of *4 . 4 so 
is either an eigenvalue or in the resolvent set. Let A > 0 be small 
enough that (14.21) holds, then any eigenfunction u of 0 satisfies all the 
assumptions of the claim above, and hence ||lnM|| < C In«4 . 4 m = 0 


so that u = 0 inside il. Clearly this implies that u = 0 in R 3 which is 
a contradiction. In the same way we deduce a uniform bound 1 /C for 
the operator V(A. 4 )~ 1 V for such small A > 0. 


Step 1.3. l! j) l 'P is well defined and bounded. Fi¬ 

nally we give a meaning to 'P(«4. 4 ) _ 1 'P, (which is required as «4 . 4 is not 
invertible on the whole space). We define it to be the strong operator 
limit of T > (A. a )~ 1 V as A —> 0. Indeed, suppose that / is fixed with 
support in Q and X n —> 0. Then we wish to compute the limit of 
for u„ = (^t 4 ") _1 7 7 / as n —> oo and show that it is independent of 
the sequence A„ —► 0. Indeed u n will satisfy 

^4-4 — X n u n T /C m n — f • 

By the same argument as before we can extract a subsequence and 
limit u £ L 6 (R 3 ;R 4 ) with convergences as in Step 1.1. In particular 
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7 7 u n —> T’u. Furthermore, any convergent subsequence u nk for any 
sequence \ n —> 0 has the same limit of Vu nk as otherwise their dif¬ 
ference would be an T 6 -eigenfunction in the kernel, which we assumed 
to not exist. We have shown that for any sequences A„ —> 0 and u n 
we can extract a subsequence converging to a limit that does not de¬ 
pend on the subsequence or the original sequence, therefore the whole 
sequence converges to this limit, which completes the proof of the ex¬ 
istence of the strong operator limit. The bound in operator norm is 
carried over by strong continuity. 

Step 2. neg(*4 4 ) = neg(«4 4 ) for all A 6 [0, A*]. A 4 is norm 
resolvent continuous in A > 0, so the only way the number of negative 
eigenvalues could change is for an eigenvalue to be absorbed into the 
essential spectrum at 0 as A —> 0. Assume this happens, then we have 
a sequence A n —> 0, a sequence of negative eigenvalues a n —> 0 and 
eigenfunctions u n which satisfy 

T A n u n T 7C u n — (7 n u n . 

By the same argument as in the previous steps, we may take subse¬ 
quences and obtain a contradiction. 

Step 3. neg(A 2 + (B 4 )*(A 4 )~ 1 B 4 ) > neg (^2 + (®°)*(«4 4 ) _1 b\) 
for all A € [0, A*]. This may be proved in the same way as Step 3 
of ILcmma 4.11 

Step 4. neg (^2 + > neg(^4 4 ). As * 8 ° = 0 and 

= 0 we have, 

neg(A2 + 

= neg(^ + 0 B?)*(.A?)- 1 ^) > neg(A?) + neg(A°) = neg(A 4 ) 
where the inequality is obtained from the assumption of the lemma. □ 

4.2.2 Finding a non-trivial kernel 

The next few steps of the proof follow those of the 1.5d case, hence we 
only provide a short overview. 

Truncation. As the domain is unbounded, each Laplacian appear¬ 
ing in the problem contributes an essential spectrum on [ 0 , oo). We 
therefore introduce a smooth positive potential function W : ffi 3 — > 
R satisfying W{x) —► oo as |x| —> oo and denote by W®" the n- 
dimensional vector-valued function with n copies of W. Then we de- 
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— A 
Mr = 


^4n 

^2,e 

~ A 

h 4 


(84V 


-a; 


-4,eJ 


'-A + eW® 3 0 

0 A - eW® 4 


A 2 

0 


0 

-A 2 


-J 


As above we can naturally define finite-dimensional operators A4 e n , 
for which we can easily prove: 

Lemma 4.7. Fixe > 0,n 6 N. Suppose that there exist 0 < A» < A* < 
— A, —A* 

00 such that neg(A4 e n ) ^ neg(A4 £ n ). Then there is a A £jn £ (A*, A*) 
for which ker(A4 e n ) is non-trivial. 

The spectrum for large A. This is again similar to the 1.5 d case, in 
particular due to the appearance of the A 2 terms. We have: 

Lemma 4.8. There is A* > 0 such that for all A > A*, e > 0, and 

—A 

n £ N, the truncated operator A4 £ n has spectrum composed of exactly 
n positive and n negative eigenvalues. In particular neg(A4 e n ) = n. 
The spectrum for small A. Again this is similar to the 1.5d case. 

Lemma 4.9. Assume that (11.101) holds and that zero is neither an 

eigenvalue of A® nor is it an L 6 -eigenvalue of A®. Then there exist 

A*,£* > 0 such that for all e £ (0,£*) there is N > 0 such that for all 
—A» 

n > N the operator A4 e n satisfies 

neg(A4 £ f n ) > neg ^ (A?) +n- neg -neg . 


5 Proofs of the main theorems 

In this section we complete the proofs of ITheorem l.ll and ITheorem 1.21 
In both settings - the 1.5d and the cylindrically symmetric - we first 
show that the results of lsection 41 imply that there exists some A > 0 
such that the equivalent problems (13.51) and (13.141) have a non-trivial 
solution (the A need not be the same in both cases, of course). Then 
we show that these non-trivial solutions lead to genuine non-trivial 
solutions of the linearised RVM in either case. 
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5.1 The 1.5 d case 


5.1.1 Existence of a non-trivial kernel of the equivalent 
problem 

Bv ILemma 4.31 and ILcmma 4.51 we have 0 < A* < A* < oo and e» > 0 
such that for any e < e* there is an N e such that for n > N e we have, 

neg(M A *J > neg(„4° + (B 0 )*^?)- 1 # 0 ) + n - neg(Al?) 

>n = neg (.M A *J, 

where the strict inequality is due to the assumption (11.91) . Fix e G 
(0, e*). By ILcmma 4.21 for each n > N e there exists A Ejn € (A*, A*) 
such that 0 G sp(At^ n ). By compactness of the interval [A*, A*] we 
may pass to a subsequence where \ e ,n k — > A E as k —> oo, for some A E G 
[A*, A*]. Bv IThcorcm 2.11 the spectra of the approximations M.^n k k 
converge in Hausdorff distance in (— K, A 2 ) to the spectrum of Al E e , 
where K > 0 is the spectral gap of the Neumann Laplacian d 2 on 
Lg(f2). This implies that 0 G sp(A / l E e ). We now repeat this argument 
to send e f 0, obtaining A G [A*, A*] with 0 G sp(.M A ). Finally, the 
discreteness of the spectrum of AT X in (—oo, A 2 ), (ILcmma 6.41) . ensures 
that 0 is an eigenvalue of AT\ i.e. A 1 a has a non-trivial kernel. 

5.1.2 Existence of a growing mode 

Now that we know that there exist some A G (0, oo) and some u = 
T 

[ip (p\ G H 2 { R) x Hg n (f2) that solve (13.51) we show that a genuine 
growing mode as defined in m really exists. To this end, we use cp, ip 
and A to define 


Ei = -d x (p E 2 = -A ip B = d x ip 

(which lie in U 1 (fl), 7J 2 (R) and il 1 (]R), respectively) and to define 
/ ± (x, v) as in ( 13 . 21 ) : 

/* = ± Vpi’ ± A(A + V ± )~ l [^{-(p + v 2 ip)\. 

Observe that / ± are both in L 2 (KxR 2 ) since and /zjf are continuous 
functions that are compactly supported in the spatial variable which 
satisfy the integrability condition ( 11 . 61 . In fact, are in the domains 
of T>±, respectively, since er- and p± are constant along trajectories 
and <p and ip are twice differentiable. 

Lemma 5.1. The functions solve the linearised Vlasov equations 

Proof. This is almost a tautology: applying the operators A + T>± 
to the expressions for f±, respectively, one is left precisely with the 
expressions (EED- □ 
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Lemma 5.2. The functions f ± belong to L 1 (R x R 2 ). 

Proof. Dropping the ± for brevity, the first term making up / is esti¬ 
mated as follows 

II Me^ll L 1 (R 3 ) llA t e||_L 2 (R 3 ) ||^’||_L 2 (R) || Me 11 (R3 ) IIMe llii(R3) ||<^|| L 2 (R) < 00 ■ 

The other terms are estimated similarly (for the terms involving the 
averaging operator this may be seen by writing the ergodic average 
explicitly isee IRcmark 3.11) or by using boundedness of the averaging 
operator on £±). This implies that f± £ L 1 (R x R 2 ). □ 

We now define the charge and current densities p and ji by 

P=J(f + ~f~)dv ji = J i>i{f + - /") dv, i = 1,2. 

Integrating in the momentum variable v alone, we obtain that 
p £ L 1 (R) as well as ji £ L 1 (R) since |{)j| < 1. In particular p,ji are 
distributions on R. 

Lemma 5.3. The continuity equation Xp + d x ji = 0 holds in the sense 
of distributions. 

Proof. This follows from integrating the linearised Vlasov equations in 
the momentum variable. Indeed, we informally have 

J (A + T)±)f ± dv = ± J [(A + D±)(pf(j> + Ppip) + X pf{-(/) + v 2 if)\ dv 

= ± J \i/j (p± + pffa) dv ± J T>± ( pf(f> + Ppif) dv = 0, 

where the Hrst term on the right hand side vanishes due to the identity 
(11.2511 and the second term vanishes since p± are even in Vi, whereas 
T>± = V\d x when applied to functions of x alone (recall that p± are 
constant along trajectories of T>±, as are pf and p±). We obtain the 
continuity equation by subtracting the ” expression above from the 
“+” expression. Owing to the low regularity of / , p and j\ this is 
true in a weak sense. □ 

Lemma 5.4. Maxwell’s equations (11.2211 hold. 

Proof. Equations (U.22bll and (U.22cl) hold due to (13.511 and the defini¬ 
tions of the operators (13.811 . Indeed, from the second line of (13.511 . we 
have 

o = ^2 Pp dv^J + J J^ptQi [h'f] dv + a 2 d> - j - 1 )(fdv 

= dl<t>+ I + pfQ± [v2^]~ pt(Q±-W) dv 
J ± 

^ dl<t> + J ( f + - f ~) dv. 
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which is (11.22el) . Similarly, (11.22bl) is obtained from the first line of 

We therefore just need to show that (11.22al) holds. However this is 
a simple consequence of (11.22cl) and the continuity equation. Indeed, 
we may first write 

—\o x Ei = A d x (f> = -A p = d x ji 

which is the derivative of (ll.22al) . Next, as (f> £ 77 2 0 (fl), its derivative 
E\ vanishes on dfl, and j\ also vanishes there due to the compact 
support of the equilibrium in 12. Thus, —A£q and j\ have the same 
derivative inside 1 1 and the same values on <912, which means they must 
be equal. □ 

Lemma 5.5. The charge and current densities p, j\ and ji are ele¬ 
ments in L 1 (ffi.) n L 2 (R). 

Proof. This follows from Maxwell’s equations and the regularity of if, 
which are in iJ 2 (R) and ff 0 2 n (Sl) respectively. □ 

This concludes the proof of IThcorcm 1.11 

5.2 The cylindrically symmetric case 

5.2.1 Existence of a non-trivial kernel of the equivalent 
problem 

The proof of the existence of a non-trivial kernel in the cylindrically 
symmetric case is in complete analogy to the one in the 1.5d case 
presented in Isubsubsection 5.1.11 and is therefore omitted. 

5.2.2 Existence of a growing mode 

Let A > 0 and u = [Ag ip A rz ] T S t/|(R 3 ;R 3 ) x H^ yl ( R 3 ) x 

H 2 Z (K 3 ;R 3 ) be such that (13.141) is satisfied, i.e. At u = 0. Let 
H 2 y ;(R 3 ; R 3 ) 9 A = Aj + A rz as in (11.321) and define 

E = —B = V x A 

(which each lie in H^^R 3 ; R 3 ) C il 1 (R 3 ;R 3 )). Furthermore, define 

/ ± = ± rpp (A • e e ) ± pfX(X + T>±y 1 (-y + A • v). 

As in the 1.5d case we begin by establishing that / ± are integrable 
and satisfy the linearised Vlasov and continuity equations. The proof 
of this result is analogous to the corresponding results in the 1.5d case, 
so is omitted. 
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Lemma 5 . 6 . The functions f± solve the linearised Vlasov equations 
(ll.29al) in the sense of distributions, and belong to L 1 (R 3 x R 3 ). Fur¬ 
thermore, the charge and current densities p and j defined by 

P=J(f + -f~)dv j= J v(f + - f~) dv, 

belong to L 1 (R 3 ) and L 1 (R 3 ;R 3 ), respectively, and satisfy the continu¬ 
ity equation Xp + V • j = 0 in the sense of distributions. 

Next we recover Maxwell’s equations from (13.141) and the continuity 
equation. 

Lemma 5 . 7 . Both the Lorenz gauge condition X ip + V • A = 0 (see 
(11.311 ) ) and Maxwell’s equations (11.261) are satisfied. 

Proof. In the same way as the 1.5d case, (ll.26al) is obtained from the 
second line of (13.141) . (13.151) and the definitions, (13.1611 . of the operators 
therein. Similarly, (ll.26bl) is obtained from the first and third lines of 

(Emem 

It remains to show that the Lorenz gauge condition holds. Using 
(11.26al) and (ll.26bl) in the continuity equation, we have, in the sense 
of distributions, 

0 = A(—A + A 2 )<p + V • [(—A + A 2 )A] 

= (—A + A 2 )[A<£ + V • A], 

As —A + A 2 is invertible, this implies that Xp + V ■ A = 0. □ 

This concludes the proof of IThcorcm 1.21 


6 Properties of the operators 

Here we gather all important properties of the operators defined in 
Isection 31 as well as the operators defined in (11.241) and (11.331) . 

6.1 The 1.5 d case 

As the only dependence on A is through the operators Q± we start 
with them: 

Lemma 6.1. In the respective spaces £± , Q± satisfy: 

(a) ||2±Hs(£ ± ) = 1- 

(b) Q± can be extended from A > 0 to Re A > 0 as a holomorphic 
operator valued function. In particular it is continuous for A > 0 
in operator norm topology. 

(c) 4sl 3 A - > oo, Q± A 1, and for u € 2)(X>±), || (Q± — l) w ||,e ± — 
\\V ± u\\ s± /A. 
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(d) As A -► 0, Q x 


± converges strongly to the projection operator Q 
defined in Definition l.f\ 

(e) For any A > 0, Q± is symmetric. 


o 

± 


Proof. ||Q^|| < 1 follows from ||(X>± + A) 1 || £± < ■pyy as zD ± is 

self-adjoint and the nearest point of the spectrum of T>± is 0. That 
ll^±ll®(£±) = 1 i s P rov ed by observing that Q±1 = 1. Part (b) follows 
from the analyticity of resolvents as functions of A. For (c) we compute 
using functional calculus for u £ l D(T>±): 


IQ+u-u L = 


£± 


/ A 2 




W-n- 1 ) 

u 

£± 

V± 



1 

A + V± 

<B(£±) 

A 

-v± 


V\ 


A 2 - Pi 


£± 


< 


< 1 • - • ||Z>±u|| 


l|Z>±«lli 


£± 


0 as A 


®(£±) 
oo 


and deduce the strong convergence Q± A 1 by the density of S5(2?±) 
in £±. 

For (d) we introduce the spectral measure of the selfadjoint operator 
—iV±, which we denote by M±. The projection onto ker(P±) is then 
= M±({0}) = f R x( a ) dM±(a) where y(0) = 1 and x( a ) = 0 when 
a ^ 0. Recall that A(A + V±)~ 1 = f R x x ia dM±(a). We compute for 
u £ £±, 


X(X + D±)- 1 u-M ± ({0})u\\ Z &± 


t——. -X(«) dM±(a)u 

A + la 




A 


A + ia 


-— x(«) 


d\\M±(a)u\\ s± , 


the last equality being due to orthogonality of spectral projections. 
This now tends to 0 as A —► 0 by the dominated convergence theorem. 
Replacing D± with —T>±, which has the same kernel, we deduce that 
A(A-Pi)- 1 A Q°±. Finally we have Q± = A(A-P±)- 1 A(A+P±)- 1 A 
(Qj_) 2 = Q± by the composition of strong operator convergence. To 
show (e) for A > 0 we simply note that D± are selfadjoint, and extend 
to A = 0 by the strong operator convergence. □ 


These results carry through to the other operators. 

Lemma 6.2. The operators jj x and B x have the properties: 

(a) For all A £ [0, oo), B x maps L 2 (R) into L 2 (f2) and J x maps 

L 2 (R) x L 2 (n) L 2 (R) x L 2 (D). 

(b) The familes {j7’ A }^ e ro i0 o) and {£> a }a<e[o,oo) are both uniformly 
bounded in operator norm. 
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(c) Both (0,oo) 9 A i —> fj x and (0,oo) 9 A 1-9 B x are continuous in 
the operator norm topology. 

(d) As A -9 0, J x - 9 . J° and B x -9 £° in f/ie strong operator 
topology. 

(e) For any A > 0 the operator J rX is symmetric. 

(f) Let V he the multiplication operator acting in L 2 (R) x Lq(A2) 
defined by 

v=\ ln °' 

L 0 In. 

where 1q is the indicator function of the set Cl. Then jj x = 

J X V. 

Proof. Part (a) is easily verifiable. We note that due to the relation 


B x = - [0 1] J x 


it is sufficient to prove the results for J x . We observe that due to the 
decay assumptions (11.61) on p±, the moment 


E 


' p± l -±^L dv 


{vY 


is bounded in L°°(R) and is real valued, so it is a bounded symmetric 
multiplication operator from L 2 (R) to L 2 (R). Next we decompose the 
second part of J x as 


E / £T±{Qi - 1)7± 
± J 




dv 


( 6 . 1 ) 


where 7± : £± x £± 9 > £± is multiplication by the vector [u 2 — l], 

and we have used the natural (and bounded) inclusions from L 2 (R) 
and Lj}(f2) into £±- Clearly 7± are bounded and we know that Q± 
have bound 1 bv lLcmma 6.11 Finally, we note that due to the decay 
assumptions on p^ and its compact support in x, that multiplication by 
pf followed by integration dv is bounded from £± to L 2 (R) and L 2 (Cl). 
Therefore jj x has a uniform bound in operator norm. Parts (c) and 
(d) then follow from the corresponding results for Q± in ILcnnna 6.11 
using (16.11) . (e) is clear from the symmetry of Q± and (16.11) . Finally (f) 
follows from the compact spatial support of p ± ,pf,Pp inside Cl. □ 

Lemma 6.3 (Properties of A x and A%)- Let 0 < A < 00 . 

(a) The operator A x is self-adjoint on Lq(CI) and the operator A\ 
is self-adjoint on L 2 (R) with the respective domains Hq n (Cl) and 
H 2 (R). 
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(b) Both [0,oo) 9 A H > Ai and [0,oo) 9 A 4 A\ are continuous in 
the norm resolvent topology. 

(c) The spectrum of A) is purely discrete. The spectrum of A^f is 
discrete and made up of finitely many eigenvalues in (—oo,A 2 ) 
and continuous (possibly with embedded eigenvalues) in [A 2 ,oo). 

(d) There exist constants 7 > 0 and A > 0 such that for all A > A, 

A ( >7, * = 1,2. 

Proof. Clearly —d 2 is symmetric. The perturbative terms are symmet¬ 
ric as well since Q± are symmetric, see lLcmma 6.11 Self-adjointness is 
guaranteed by standard arguments, such as the Kato-Rellich theorem. 

Let us prove (b), considering first Af). It is sufficient to prove that 
(A2 — *) _1 —t (AS — *) _1 in operator norm as A — > a (with A,ct > 0). 

We use the second resolvent identity to obtain 

- i )- 1 - (AS - i r 1 = (A - i)-\AZ - A$)(A% - i )- 1 

= (-42 - - A 2 ) - (J& - JA))(AZ - i)- 1 

where J{\ is the upper left component of written in block matrix 
form. Hence, as the resolvents are each bounded in operator norm by 

1 , 

11 0^2 — *) ~ (A 2 ~ ll<B(L 2 (R)— I 17 ^A l + ll (3\\ ~ Jll)(A 2 — *) ||<g(i2 

It thus suffices using iLninma 6.2l fl to show that (3°\ — 3\i)V(A^ — 
i ) _1 —> 0 in operator norm, where V is the multiplication operator 
on L 2 (R) given by the indicator function of the set Q. V is relatively 
compact with respect to — d 2 by the Rellich theorem, and hence also 
relatively compact with respect to A% as it also has the domain H 2 (R.) 
by part (a). Hence V(A% — *) _1 is compact, which allows us to upgrade 
the strong convergence jf\ A given by ILcnnna Ol to operator 
norm convergence. The norm resolvent continuity of A 2 follows. The 
proof for Ai is analogous, but lacking the |cr 2 — A 2 | term. 

Part (c) is simple: both operators have a differential part (Lapla- 
cian) and a relatively compact perturbation. Hence both conclusions 
follow from Weyl’s theorem [5J IV, Theorem 5.35]. The finiteness of 
the discrete spectrum below the essential part in the case of A 2 is a 
consequence of ILcnnna 2.21 For part (d), we begin with A). Fix an 
arbitrary h £ Hq u (TI), then we compute, 

(Aih,h) L 2 {n) = \\dxh\\ 2 L 2 {n) -J2ff hfj,f(l- Q±)hdvdx 
Now we note as h £ H 2 n (Q), h is in 33(X>±) when interpreted in £±. 
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We now use lLcmma B.lf cl to estimate, 

(A*h, h) LUn) > \\d x h\\l i{ n) \H/ W± IL~ ( n X R3) ll^ll £± M s± 




> K\ 


Hig(fi) 


1 - 


C 


Vk A 


where K is the spectral gap of the Laplacian on the bounded domain 
SI, and we have used 

\[D±h\\ £ ± = Jj w ± \v 1 d x h\ 2 dvdx < \\d x h\\ 2 L 2 ( n) sup J w ± \vi\ 2 dv 

and the natural bounded inclusions from Lq(£1) into £±. We now just 
take A > C/VK. 

For A% the proof is easier due to the A 2 term. For h £ iJ 2 (R) we 
compute, using the formulation (13. Sail . 


(A$h,h) LHm = \\d x h\\l*<n + ^\\hf LHm - 
> (A 2 - C') I 


h 


'h 

0 

? 

0 


i 2 (R)xig(H) 


li 2 (R) 


where we have used that the uniform bound in operator norm of j 
given bv ILcmma 6.21 We now take A > \[C''. □ 

Lemma 6.4 (Properties of A4 a ). For each A £ [0, oo), the operator 
A4 a is self-adjoint on L 2 (R) x Lq(S7) with domain H 2 ( R) x Hq n (S7). 
For any A > 0, the operator has essential spectrum [A 2 ,oo). The 
family {A4 A }^ e [o j0 o) is continuous in the norm resolvent topology. 

Proof. The proof essentially mimics (and uses) the preceding proofs, 
and is therefore left for the reader. □ 


6.2 The cylindrically symmetric case 

As many of the proofs are the same as in the 1.5 dimensional case 
above, we give the details only where they differ. 

Lemma 6.5. In the respective spaces 9t±, Q± tSym and Q± skew satisfy: 


(a) 


Qi 

~^'±,sym 


®(W±) 

‘A 


= 1 and 


O x 

,skew 


< 


<B(OT±) 


( b ) Q±,sym and Q±,skew can be extended from A > 0 to Re A > 0 
as holomorphic operator valued functions. In particular they are 
continuous for A > 0 in operator norm topology. 

(c) As R 9 A — )■ oo, Q± ySym —> 1, and for u £ S(2?±), we have the 


bound 


( Q±,sym - !)* 


< 




V±U 




/A. 
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(d) islgA 
bound 


oo, Q 


±,skew 


0, and for u £ ®('D±) ; we have the 




< 




V±u 




/A. 


(e) As A —> 0, Q± aym converge strongly to the projection operators 
Q± defined in \Definition 1.5\ 

(f) As A —> 0, Q± skew converge strongly to 0. 


(g) For any A > 0, Q± tSym are symmetric and Q± skew are skew- 
symmetric. 


Proof. The claims about Q± sym may be proved in the same way as 
those in lLcmma 6.11 For (a), the spectral theorem applied to the self- 
adjoint operators — iD± implies that Q± skew are unitarily equivalent 
to a multiplication operator, so that 


Q 


±,skew 


s(ai±) 


—ia A 
A 2 + a 2 


< 

L“(sp(if5 ± )) 


—ia A 
A 2 + a 2 


1 

2 ' 


The proof of (b) follows, as in the proof of lLcmma 6.11 from the holo- 
morphicity of the resolvent. For (d), we let u £ 3?(2?±), and then for 
A > 0 we have, 


Q 


A , 
±,skew L 


n± 


< A 


1 

< - 

- A 


(A 2 + 
V±u 


vlr 1 


<n ± 


T>±u 

<s(m ± ) 

0 as A —> oo. 




The strong convergence to 0 then follows from the density of T>(fD±). 

For (f), we repeat the proof of lLcmma 6.11 noting that it is shown that 
A(A + V ±)~ 1 -> Q° ± and Q± iSym -> Q± as A -*■ 0. That Q± tSkew —> 0 
as A —> 0 now follows from the identity, valid for all A > 0, 

A(A + X>±) 1 = Q±, sym + Q± tS keW 

Finally, (g) is a consequence of lLcmma 3.21 □ 

~ A ~ 

Lemma 6.6. The operators jj , and Bf for i = 1,2, 3,4 have the 
properties: 

(a) For all A £ [0,oo), Bf £ 03 (l 2 (R 3 ; K 3 ), L 2 yi (R 3 )) , B\ £ 03 (L^(R 3 ; R 3 ), L 2 2 (R 3 ; R 3 )), 

£ B (.L^(R 3 ),L 2 2 (R 3 ;R 3 )) andB, £ 23 (l 2 (R 3 ; R 3 ), L^(R 3 ) x L 2 2 (R 3 ;R 3 )) 
with bounds uniform in A. 

(b) Each of (0, oo) 9 A J and (0, oo) 9A4 Bf, i = 1,2, 3,4 are 
continuous in the operator norm topology. 

(c) As A —^ 0, jj —> J , B 4 — > B\, 0 and B$ 0 in the 

strong topology. 

(d) For any A > 0 the operator J is symmetric. 
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(e) Let "P be the multiplication operator acting in L 2 cyl 
L 2 2 (R 3 ;R 3 ) defined by 


V = 


In 

0 

0 


0 

In 

0 


0 ' 

0 

In. 


x 


where 1q is the indicator function of the set LI. 


J V. 


Then J = 


Proof. That the operators map the corresponding spaces to each other 
may be verified directly from (13.161) . noting in particular the notation 
vg = egi>g and v rz = e r v r + e z v z . As in the proof of ILcmma 6.21 the 
uniform (in A) bound on the operator norms may be obtained using 
the decay assumptions on the equilibrium and the uniform bound on 
the norms of Q± sym and Q± skew given by ILcmma 6.51 In the same 
way (c) and (d) follow from the corresponding results in ILcmma 6.51 

~ A 

To show the symmetry of J for A > 0 we use the block matrix 
form, noting that Q± tSym appears on the diagonal, and that the off 
diagonal entries have Bf and their adjoints in the appropriate config¬ 
uration. Then we extend to A = 0 by the strong convergence. As 
in ILcmma 6.21 (e) follows from the spatial support properties of the 
equilibrium. □ 

Lemma 6.7 (Properties of A) , A%, A% and A 4 ). Let 0 < A < oo. 

(a) The operator A) is self-adjoint on L 2 yl ( R 3 ), the operator A\ 
is self-adjoint on Lg(R 3 ;R 3 ), A$ is self-adjoint on L 2 2 (R 3 ;R 3 ) 
and A 4 is self-adjoint on L 2 yl (R 3 ) x L^ 2 (R 3 ;R 3 ) with the respec¬ 
tive domains H 2 yl ( R 3 ), Hg{ R 3 ;R 3 ), U^ 2 (R 3 ;R 3 ) and H 2 yl { R 3 ) x 
LT r 2 2 (R 3 ;R 3 ). 

(b) The mappings [0, oo) 9 A A 4 , [0, oo) 9 A i —> Afj, [0, oo) 9 A9 
A 3 and [0, 00 ) 9 A 1-9 A 4 are continuous in the norm resolvent 
topology. 

(c) The spectra of A 4 , A%, A 4 and A 4 are discrete and finite in 
(— 00 , A 2 ) and continuous (possibly with embedded eigenvalues) 
in [A 2 , 00 ). 

(d) There exist constants 7 > 0 and A > 0 such that for all A > A, 

A( > 7 , * = 1,2,3,4. 

Proof. The proof for each of A(, i = 1 , 2, 3,4 is analogous to that of 
ILcmma 6.31 for A%- We omit the details. □ 

Lemma 6.8 (Properties of A4 ). For each A € [0, oo), the operator 
A4 is self-adjoint on Lg(R 3 ;R 3 )xL 2 yl (R 3 )xL 2 z (R 3 ;R 3 ) with domain 
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7/|(R 3 ; R 3 ) x H% yl ( R 3 ) x 7J 2 Z (R 3 ; R 3 ). For any A > 0, the operator Ai 

has essential spectrum (—oo, — A 2 ] U [A 2 , oo). The family {Ai }Ae[o,oo) 
is continuous in the norm resolvent topology. 

Proof. This is again analogous to the previous proofs, and is left to 
the reader. □ 


7 Construction of example equilibria 

In this section we will provide (non-optimal) existence criteria for com¬ 
pactly supported equilibria of the 1.5d system which can be written in 
the form (11.16[l - (11.17fl . Examples of equilibria in the cylindrically sym¬ 
metric case were already provided in [9]. 

Proposition 7.1 (Existence of confined equilibria). Let R > 0, a > 2 
and C R 2 be bounded domains. Then there are constants c, C > 0 
such that if two functions p ± {e ± ,p ± ) € Co(R 2 ) with support in A ± 
satisfy 

I < c(l + |e ± |)“ a 

and a function ip ext € U 2 oc (R) satisfies 

\xjj ext (x)\ > C{1 + \x\ 2 ) for \x\ > R 

then there are potentials <j>°,if° € H? oc ( R) such that {p ± {e ± ,p ± ), </>°, if 0 , ip ext , (f ext 
0) is an equilibrium of the 1.5 d relativistic Vlasov-Maxwell equations 
(11.151) with spatial support in [~R,R], where the relationship between 
(x,vi,v 2 ) and (e^p 31 ) is as defined in (11.171) . 

Remark 7.1 (Trivial solutions). Of course, the result does not say that 
the obtained equilibrium is not everywhere zero. This may be ruled out 
by choosing p± and ip ext in such a way that (for example) p ± {x = 

0,u = 0) > 0 if<p°,il> 0 = 0. Let us sketch the argument. Recall that we 
write f 0i± (x,v) = p ± {{v) ±(f 0 (x),V 2 ±ip 0 (x)±ip ext (x)) = p ± {e ± ,p ± ). 

If / 0,± (a;,n) = 0 for all (x,v) then p,ji = 0 and (f 0 ,^ 0 = 0 for all x. 

Therefore e± = (v) and p± = v 2 ± ip ext {x), and 

/ o,± (0,0) = p ± (l, ±i/> ext (0))■ 

The RHS is something we can ensure is positive by choosing A , pr 1 
and i/j ext appropriately. Under this appropriate choice one obtains a 
contradiction. 

Proof of \Proposition 7~1\ Given two elements p , j 2 € L 2 (R) with com¬ 
pact support, we define 

(j)° = G * p, ip° = G*j 2 , 
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where G(x) = — \x\/2 is the fundamental solution of the Laplacian in 
one dimension. [We note that one expects j\ to vanish for an equi¬ 
librium, due to parity in iq, hence it does not appear in the setup.] 
Thus we define e ± = e ± (x, vi,v 2 ) and p± = p ± (x, v±, V 2 ) via the usual 
relations (11.171) . which we recall for the reader’s convenience: 

e ± (x, v ) = ( v} ± (/)°(x), P ± (x, v) = V 2 ± tp°{x) ± ip ext (x) 

(c jf xt is zero). We let T : L 2 (R) 2 —> L 2 (M) 2 be the (non-linear) map 
defined by 


T 


' p 

_ f 

'1' 

_h_ 

-J 

y 2 


(d {e ,P + ) ~ P (e ,p ))dv. 


A fixed point of J- is the charge and current densities of an equilibrium 
(p ± (e ± ,p ± ), </>°, ?/>°, ip ext , (j) ext = 0). We define X C L 2 (R) 2 as 


X = {(p,j 2 )£L 2 


: both supported in [— R, R] and bounded by C'} 


for a positive constant C' to be chosen. This set is clearly convex. We 
will show that for c > 0 sufficiently small and C > 0 sufficiently large, 
.F is a compact continuous map X ^ X and thus, by the Schauder 
fixed point theorem, has a fixed point. 


Step 1: Compact support. C' and C can be chosen so that T 
maps X into functions supported in [— R, R], 

For ( p,j 2 ) G X and |x| > R , we have, 

|/(x)| = \(G*p)(x)\ <C' f \G(x—y)\ dy = [ \x-y\dy = C'R\x\ 

J-R £ J-R 

and the same bound holds for ijr. This allows us to control v 2 using 
(X and x. Indeed, 

|-U 2 1 < (u) = e ± =F < |e ± | + |</>°(a;)| < |e ± | + C'R\x\. 

Which gives the following lower bound for |p ± | + |e ± | in terms of x: 

IP *I + \e*\ = \V 2 ± 1>\x) ± r xt (x)\ + ^ > r xt (x) - M - \^°(x)\ + |e ± | 

> ip ext {x) — |e ± | — 2C"i?|a;| + |e ± | 

> C(1 + |x| 2 ) — 2C , i?|a:|. 

By taking C' small enough and C large enough we can ensure that 
if |x| > R then (e^^p^) lie outside any disc in K 2 , and in particular 
outside A , where are supported. This proves the claim. 
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Step 2: Uniform L°° bound. C' and c can be chosen so that 
T maps to functions with L°° norm smaller than C'. 

Estimating </>°( x) for \x\ < R 

|</>°0)| < y j Jx-y\dy = -^-(x 2 + R 2 ), 


we take C' small enough (recall it was already taken to be small in 
the previous step, hence we may require it to be even smaller) so that 
|</>°(x)| < 3/4 for |x| < R. Now the decay assumption on p ± allows us 
to show a uniform bound on |J r i(/9, j 2 )(x)| in |x| < R: 


1 \P ± ( e± ^P ± )\ dv < E 1 + C \ e ±\)oc dv 


^E 


(1 + {v) - |0°(x)|)‘ 
2 c 


■ dv <Y, 


(! + <«)-!)' 


■ dv 


(l + W)' 


■ dv = C c < oo. 


(7.1) 

We can bound \F 2 (p,j 2 )(x)\ in the same way as |0| < 1. Finally we 
choose c so that C"c < C'. 


Step 3: Uniform L°° bound on the derivative. There is a 
constant C'" such that for any (p, j 2 ) € X, we have \\d x F\{p, j 2 )\\ LOO ,_ R R ] < 
C'" and \\d x F 2 (p,j 2 )\\ Lao[ _ RtR] < C'". 

We compute for F\, and note that F 2 is analogous. Using the chain 
rule, we have 

dx J (p + {e + ,p + ) - p~{e~,p~)) dv = 

{d x 4>°) J(p,t(e + ,p + ) + p~(e~,p-))dv 

+ (d x tp° + d x i/j ext ) J(p+(e + ,p + ) +p~{e~,p~))dv. 

The two integrals are bounded uniformly in x by the arguments in Step 
2 using the corresponding assumed bounds on p,f and p± respectively. 

As the external field (f> ext lies in Hf ( R), its derivative d x (j> ext lies in 
i7 1 ([— R, i?]) and is bounded in L°°([—R, i?]) by Morrey’s inequality. It 
remains to bound d x (jr and d x ip° uniformly for all x € [— R,R]. These 
are controlled directly using the Green’s function G{x ) and uniform 
bounds of Step 2. Indeed, 

\{d x 4>°){x)\ = K^G) * p)(x)| < y / I sign^ ^ y )I dy < C'R 
and the computation for d x ip° is identical. 
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Step 4: T is a compact continuous map from X to X. 

Steps 1 and 2 imply that X(X) C X. Step 3 and the Rellich theo¬ 
rem imply that X(X) is relatively compact in X. It remains to show 
that T is continuous. This may be shown using dominated conver¬ 
gence and the bounds in Step 2. Indeed, suppose that (p", ) C X 

is a sequence converging to (p, J 2 ) € X strongly in L 2 (R) 2 . We shall 
show that /i(p n ,j - 2 ) —> Xi(p,j 2 ) in L 2 , the result for T 2 is analo¬ 
gous. By Step 2 and dominated convergence it is enough to show 
convergence pointwise, i.e. for each x G [—R,R]. Next, by (17.11) and 
dominated convergence again, it is sufficient to show that the corre¬ 
sponding densities /i ± (e ± ,p ± ) converge pointwise in (x, v). Continuity 
of pA reduces this to showing pointwise convergence of the correspond¬ 
ing microscopic energy and momenta and pA The definitions of 
these quantities imply that it is enough to show that the correspond¬ 
ing electric and magnetic potentials 4>°’ n and tp°’ n converge pointwise. 
Finally, as the potentials are (p™,)?) convolved with G(x) = — \x\/2, 
the convergence (p n , jj) —> (p,j 2 ) in L 2 ([—R,R]) 2 gives the required 
pointwise convergence. 

This concludes the proof. □ 
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